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ABSTRACT 


We  consider  quadrature  formulae  (q.  f. )  for  the  numerical  evalua¬ 
tion  of  the  Fourier,  cosine,  sine,  and  Laplace  transformations.  Let 

S  denote  the  clas  s  of  spline  functions  of  degree  n- 1  defined  on 
n 

the  real  line  and  having  simple  knots  at  the  points  v  +  ~  for  all 
integers  v .  This  means  that  S(x)  c  Sr  provided  that  S(x)  c  C  ~ 
and  that  the  restriction  of  S(x)  to  any  interval  between  consecutive 
*  knots  is  a  polynomial  of  degree  not  exceeding  n  -  1. 

t  In  Part  I,  we  consider,  for  n  a  positive  Integer,  a  q.  f.  of  the 

form 

/  ffxle^dx  =  £  H<n)f(v)  +  Rf 

-00  v  =-  00  ' 

/.I 

where^  for  fixed  t,  the  coefficients  H'  ;  are  bounded.  We  show  that 

v>  ‘ 

among  all  such  q.  f. ,  there  is  a  unique  formula  with  the  property  of 
being  exact,  i.  e.  the  remainder  Rf  =  0,  whenever  f(x)  c  Sn  fl  L^R). 

We  exhibit  the  explicit  formula  for  arbitrary  step  length  h  and  give  a 
useful  bound  on  the  remainder  Rf  when  n  is  even. 

In  Part  II,  we  discuss  thu  cosine  and  sins  transforms,  using  deriva¬ 
tive  data  at  the  origin.  For  the  cosine  case,  we  consider  q.  f.  of  the 


form 


00  00  f2^] 

/  f(x)oos  xt  dx  *  2  H*n)  f(v)  +  £  y2j”1)(0)  +  Rf 

0  v=0  V»X  J=1 

where,  fox*  fixed  t.  again  the  coefficients  are  bounded,  or  the 

v,  t 

similar  case  when  f^(0)  (j  s  0,1,..., [^~ ])  is  known.  We  find 
that  among  all  such  q.  f.  there  is  a  unique  one  that  is  exact  when 
f(x)  c  f!  Lj(R+).  We  exhibit  the  explicit  q.  f.  for  arbitrary  n,  but 
have  a  proof  only  for  specific  n. 

In  Part  III,  the  Laplace  transform  case,  we  use  the  weight  func¬ 
tion  e~Xp  instead  of  cos  xt  or  sin  xt  as  in  Part  II,  but  otherwise 
proceed  in  much  the  same  spirit  as  Part  II.  Part  IV  contains  expressions 
for  the  remainder  or  error  for  the  q.  f.  in  the  first  three  parts  and  explicit 
error  bounds  for  the  approximations  of  the  first  tv/o  parts.  Two  compu¬ 
tational  examples  are  also  included. 

We  actually  use  three  different  approaches  to  construct  our  q.  f. : 
we  either  integrate  an  appropriate  spline  interpolant  to  f(x),  require 
our  q.  f.  to  be  exact  for  a  particular  sequence  of  so-called  B- splines, 

f 

or  utilize  a  particular  monospline.  In  any  case,  the  generality  and 
utility  we  achieve  is  due  to  the  form  of  the  splines  we  use,  in  particu¬ 
lar  to  the  components  of  these  splines,  the  so-called  B- splines. 
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INTRODUCTION 


In  [10]  I.  J.  Schoenberg  generalized  the  construction  of  best  quad¬ 
rature  formulae  in  two  ways.  He  discusses  integrals  with  an  arbitrary 
pre-assigned  weight  function  opening  up  the  possibility  of  constructing 
quadrature  formula  (q.  f. )  for  the  numerical  evaluation  of  Laplace  trans¬ 
forms,  Fourier  integrals,  and  other  special  integral  transforms.  We 
pursue  this  possibility  here;  in  particular  we  wish  to  discuss  approxi¬ 
mations  to  the  integrals 


r  ixt 

(1)  /  f(x)elx  ax, 


(2)  /  f(x)cos  xt  dx, 
0 


(3)  /  f(x)sin  xt  dx, 

0 


w 

(4)  /  f(x)e”xpdx 
0 


which  are  the  Fourier,  cosine,  sine,  and  Laplace  transformations, 
respectively. 

In  1949,  A.  Sard  generalized  the  Newton-Cotes  q.  f.  as  follows:  let 

1  £  m  £  n  +  l  and  let 

n  n 

(5)  /  f(x)dx  =  J  H(m)  f(v)  +  Rf 

0  v=0  v’ n 

be  the  formula  exact,  i.e.,  Rf  =  0,  if  f(x)  e  n  ,  ,  the  class  of 
polynomials  of  degree  not  exceeding  m-1,  and  such  that  the  func¬ 
tional  Rf  when  written  in  Peano- fashion  as  an  integral  of  the  form 

/  k(x)f  (x)dx  has  the  kernel  k(x)  with  least  L-norm.  It  was 
0  2 
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shown  by  Schoenberg  [8]  that  we  can  say  the  following:  8ard's 
q.  f.  (5)  is  uniquely  characterized  by  the  requirement  of  being  exact, 
hence  Rf  =  0,  for  the  elements  of  the  class  (0, 1, .  „ . ,  n)  of 
natural  splines  of  degree  2m  -  1  having  the  knots  0, 1, . . . ,  n. 

The  term  "natural"  indicates  that  the  degree  of  the  polynomial  com¬ 
ponents  of  the  spline  function  should  drop  from  2m  -  1  to  m-1 
in  each  of  the  two  intervals  (-oo,  0)  and  (n,  oo). 

In  [10]  Schoenberg  discusses  q.  f.  of  the  form 

(6)  / w(x)f(x)dx  =  £  H(mJf(v)  +  "§l  |W  ,«>«,)  +"£  c[m^f(,,(n)  +  Rf 

0  v=0  V»n  J*1  J=1  ,,n 

where  w(x)  is  an  arbitrary  preassigned  weight  function  and  such 

that  the  q.  f.  (6)  is  exact  for  II  ,  and  has  the  property  that  the  a»~ 

m-1 

sedated  kernel  k(x}  of  the  functional  Rf  has  least  L^-norm.  This 
q,  f.  he  shows  is  uniquely  characterized  by  requiring  Rf  *  0  if  f  is 
a  spline  function  (not  natural)  of  degree  2m  -  1  having  the  knots 
1,  ...,n-l. 

In  the  paper  [12],  Schoenberg  discusses  infinite  analogues  of 

Sard's  q. f.  (5)  for  the  real  line  F  and  the  half-line  F+  or  (0, oo). 

We  first  consider  the  entire  line,  the  so-called  cardinal  case  when 

all  integers  v  are  nodes  of  the  formula.  Let  n  be  a  natural  number 

and- let  S  denote  the  class  of  spline  functions  of  degree  n  -  1,  or 
n 

order  n,  defined  on  the  real  line  and  having  simple  knots  at  the 
integers  v  if  n  is  even,  or  at  the  halfway  points  v  +  —  if  n  is 


i. 


odd.  This  ii,aans  that  S(x)  <  SR  provided  that  S(x)  c  C^2  (for 
n  *  1  this  condition  la  vacuous),  and  that  the  restriction  of  S(x):  to 

i 

any  Interval  between  oonsecutive  knots  is  Identical  with  a  polynomial 
of  degree  not  exceeding  n  -  1,.  Such  functions  are  called  cardinal 
spline  functions. 

.  i 

Lemma  1  below  (§1)  shows  that 

00 

(7)  S(x)  «  8  (1  L.  (F)  implies  that  TV  |S(v)|<oo. 

n  l 


v=-oo 


Let  n  be  even,  say  n  =  2m,  and  consider  a  q.  f.  of  the  form 


(8) 


/  f(x)dx  =  £  H*2m)  f(v)  +  Rf 


—  00 


-00 


r«m) 


where  the  numerical  coefficients  IT  K  satisfy  the  condition  that 


(9) 


lH*2m)  |  <  K  for  all  v  and  some  appropriate  K  . 


The  Implication  (7)  shows  that  the  functional  Rf  is  well-defined  by 

I 

(8)  if  f(x) «  S2m  D  LjtF).  One  of  the  results  of  [12]  is  as  foljowst 


Among  all  quadrature  formulae  (8).  (9)  the  q.  f. 


(10) 


00  *» 

/  f(x)dx  =  £  f(v)  +  Rf 


•00 


•00 


is  characterized  byjhe  requirement  that  Rf  =  0  il  f  c  S2m  n  L^F). 

Observe  that  (10)  is  none  other  than  the  Euler- Maclaurin  q.  f, 
oo  oo  rn  oo 

J  f(x)dx  =  £  /  Bm(x)f(m)(x)dx, 

-oo  -oo  -oo 


th 


where,  if  B  (x)  denotes  the  m  Bernoulli  polynomial,  we  have  de¬ 
ni  ' 

fined  B  (x)  to  be  its  periodic  extension  of  period  1  from  the  interval 
m 

[0,1]. 
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In  Part  I,  We  consider  the  analog  of  q.  f.  (6)  for  the  entire  line 

ixt 

]R  and  we  take  w(x)  =  e  ,  that  is,  we  discuss  approximations  to  the 

i  • 

t 

i 

general  Fourier  transform  (1).  Let  n  be  any  positive  integer  and 

t  «  • 

consider  a  q.f.  of  the  form 

00  00 

(11)  /  f(x)eixtdx  =  £  H(n)  f(v)  +  Rf 

V  •  v 

-00  -  00  9 

where  the:  coefficients  H^.  satisfy  the  condition  that 

v,t 

i 

(12)  |H(n).J  <  K  for  fixed  t,  for  all  v  and  some  K. 

v* 1 

i 

;  Note  that  the  coefficients  are  now  functions  of  t.  Again,  the 

result  (7)  shows  that  the  functional  Rf,  now  given  by  (11),  is  well- 

I 

defined  if  f(x)  «  sn  n  L^R).  We  have  the  following 

i 

Theorem  1.  Among  all  quadrature  formulae  of  the  form  (11).  (12), 
there  is  a  unique  formula  with  the  property  of  being  exact,  i.  e. .  Rf  =  0, 

;  whenever  f(x)  <  Sn  D  Lj(R). 

1 

This  q.  f.  (11)  could  also  be  obtained  by  using  Newton's  funda¬ 
mental  idea;  assuming  the  function  f(x)  to  be  given  numerically  at 
equidistant  points  of  step  1,  including  the  origin  0,  we  interpolate  f(x) 
by  a  function  S(x)  at  these  points,  and  then  construct  the  Fourier 
transform  of  S(x).  This  idea  has  been  used  before,  and  often,  for  the 
integrals  (1)  -  (4)  [5]. 

In  fact,  for  n  =  2,  the  case  of  linear  spline  interpolation,  the 

i 

q.f.  (11)  can  be  found  in  [5,  pp.  22,  23],  But  for  n  >  2,  similar 


-4- 
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q.  f.  (11),  (12)  have  not  previously  been  developed.  The  generality 
and  utility  we  achieve  is  due  to  the  form  of  the  interpolatin'-  'unc¬ 
tions  we  use,  i.  e. ,  the  splines,  and,  in  particular,  to  the  compo¬ 
nents  of  these  splines,  the  so-called  B- splines. 

In  the  paper  [12],  Schoenberg  also  considers  the  analog  of  Sard's 

q.  f.  (5)  for  the  half-line  F+.  Let  st  denote  the  class  of  functions 

2m 

S(x)  satisfying  the  following  four  conditions. 

1*  S(x)  c  C2m“2(R) 

2#  S(x)  c  n_  ,  in  each  interval  (v,  v+  1)  for  v  =0,1,... 
cm-i 

3*  S(x)  c  n  .in  the  interval  (-»,  0) 
m- 1 

4*  S(x)  c  L1(R+). 

We  now  consider  a  q.  f.  of  the  form 
00  00  l?.  ) 

(13)  /  f(x)dx  =  £  «*)  + Rf 

0  v=0 

whose  coefficients  satisfy  the  condition  that 

(14)  |H*2ir)|  <K  for  v  £  0  and  some  K. 

CO 

By  Lemma  5  of  [12],  S(x)  «  S*  implies  that  ^  |S(v)|  <  oo  so  that 

feUl  A 

v  =0 

the  q.  f.  (13)  is  applicable  whenever  f(x)  e  st  . 

2m 

In  [12],  Schoenberg  proved  the  following 

(lii)  Among  all  qt  f.  of  the  form  (13).  (14),  there  is  a  unique 

formula  with  the  property  of  being  exact,  i.  e. .  Rf  =  0, 
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4 

whenever  f(x)  <  S_  . 

*~“  /<m 

In  Part  II,  we  consider  the  analog  of  q.  f.  (6)  for  the  half-line 

4 

R  and  we  take 

(16)  w(x)  =  cos  xt  or  w(x)  =  sin  xi . 

We  want  a  q.  f.  of  the  form 

00  00  jiT—  1 

(17)  /  w(x)f(x)dx  =  Yj  H(2m)f(v)+  £  bJ2"0  f(J)(0)  +  Rf 

0  v=0  v,t  j=l 

where  w(x)  is  given  by  (16)  and  the  coefficients  H^2^  satisfy  the 

v»  * 

condition  that 

(18)  |h'  'i  <  K  for  fixed  t,  for  all  integers  v^O  and  some.K„ 

v» 1 

Note  that  for  m  fixed,  H^2m^  is  a  function  of  t.  Lemma  1  below 

v,t 

(§  1)  shows  that 

00 

(19)  S(x)  c  S  n  L1(R+)  implies  that  Y  IS(v)|  <  00 

v=0 

so  that  the  functional  Rf  is  well-defined  if  f(x)  e  S,  fl  L,(R+). 

2  m  1 

Similar  to  (15),  then,  for  our  endpoint  derivative  case,  where  w(x) 
is  given  by  (16),  we  have 

Theorem  2,  Among  all  quadrature  formulae  of  the  form  (17),  (18), 

there  is  a  unique  formula  with  the  property  of  being  exact,  i.  e. . 

Rf  =  0.  whenever  ffx)  «  S_  n  L.  (R+). 

1  -  2m  1 

We  also  consider  q.  f.  of  the  form 


1 


1 


* 


i 

i 

\ 

I 

1 

1 

( 


i 

5 


> 


* 
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(21) 


B(n)  f(2J) 


(0)  +  Rf 


/|*% 

where  the  coefficients  IT  ’  again  satisfy  the  condition  (18).  Lemma 

v» 1 

1  will  again  assure  us  that  the  functional  Rf  given  by  (20)  or  (21) 

is  well-defined  if  f(x)  <  S  0  L  (F+).  The  form  of  the  H(n*  and  the 

«  *  v,  t 

®  j  |  (20)  and  (21)  is  particularly  simpler  than  the  corresponding 

form  of  the  coefficients  for  the  q,  f.  (17)  and  for  this  reason  we  shall 

consider  the  q.  f.  (20)  and  (21)  first.  Our  approach  is  the  following: 

once  the  existence  and  unicity  of  the  q.  f.  have  been  established,  we 

shall  exhibit  the  q.  f.  (20)  and  (21)  that  satisfy  the  requirement  (18) 

and  show  that  they  are  exact  whenever  f(x)  c  n  L^R*)  . 

We  could  also  obtain  the  q.  f.  (20)  and  (21)  by  constructing  the 

cosine  or  sine  transform  of  the  appropriate  spline  interpolant.  Closest 

tc  this  point  of  view  is  the  paper  [lj  in  which  Einarsson  approximates 

integrals  of  the  form 


b  b 

(22)  /  f(x)  cos  wx  dx,  /  f(x)  sin  wx  dx 

a  a 

as  follows:  f(x)  is  interpolated  by  a  cubic  spline  with  equidistant 
knots,  the  interpolation  being  at  the  knot  s,  while  the  values  f'(a) 
and  l»(b)  are  matched  by  the  cubic  spline.  Then  he  takes  the  trans- 
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form  of  the  spline.  We  could  adopt  this  method  al90,  and  through  the 

use  of  B-  spline s,  achieve  greater  generality  than  Einarsson.  However, 

we  do  not  follow  this  approach  because,  in  the  general  case,  this 

method  does  not  lead  to  the  coefficients  ,  Bjn?  of  (20)  and 

v,  t  j,  t 

(21)  in  a  very  simple  form. 

In  Part  III,  we  consider  the  Laplace  transform  (4)  and  establish 
the  following 

Theorem  3.  Among  all  q.  f.  of  the  form 

(23)  /  f(x)e"Xpdx  =  £  H  m  f(v)  +  £  B*m)  f(,,(0)  +  Rf 

0  v=0  v,p  J=1  J»p 

whose  ooerfl dents  satisfy  the  condition 

(24)  jH^  |  <  fyi  V  P  fixed,  for  all  v  H  0,  and  some 
K,  some  p  >  1, 

there  is  a  unique  formula  with  the  property  of  being  exact  whenever 
f(x)  is  a  cardinal  spline  function  of  degree  2m-  1  such  that 
f(x)  =  0(x8)  as  x-  no  for  some  s  =  0. 

We  do  this  in  the  same  way  we  prove  Theorem  2,  by  using  a 
generating  function  approach  similar  to  that  used  in  [13], 

Part  IV  contains  expressions  for  the  error  for  the  approximations 
we  make  in  the  first  three  parts  and  estimates  of  error  bounds  for  the 
first  two  parts.  We  acquire  these  expressions  by  showing  that  we 
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4 


could  have  constructed  our  q,  f.  in  still  a  third  way.  This  third  ap¬ 
proach  utilizes  a  particular  mono  spline,  related  to  the  so-called 
Rodrigues  function  of  [10]. 


j 

i 

I 

*. 

i 

\ 

i 


> 


t 
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I.  THE  FOURIER  TRANSFORM 


1.  Preliminaries.  We  first  recall  some  known  definitions  and 
results  [7].  Let  n  be  a  natural  number  and 

(l.l)  MM  =  Mn(x)  =  6nx^1 

where 


f  x  if  x^O 

X+  \o  ifx<0 

and  where  6n  stands  for  the  usual  symbol  of  the  n**1  order  central  dif¬ 
ference  of  step  equal  to  1.  Mr(x)  is  a  spline  function  of  degree  n  -  1 
having  as  knots  'lie  points  v  (v  integer),  or  v  +  ~  ,  depending  on 
whether  n  -  1  is  odd  or  even.  Mn(x)  is  positive  in  the  interval 

(-  Tn,  -rn)  and  vanishes  elsewhere,  and  evidently  M  (x)  «  S  .  It 

n  n 


has  the  following  Fourier  transform: 


00 


(1.2)  /  Mn(x)e  <ix  =  «|»n(t) 
-00 

where 

(1.3) 


2  sin£ 


Mn(x)  is  called  a  central  B- spline  or  basis  splice  because  of  the  fol¬ 
lowing  property:  If  &(x)  <  S^,  then  S(x)  admits  a  unique  representa¬ 


tion  of  the  form 

(1.4)  S(x)  =  J  C  M  (x-  v) 

<-»  v  n 


oo 


-00 
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and  conversely,  any  such  series  with  arbitrarily  prescribed  {C^  } 

V 

converges  and  defines  a  cardinal  spline  function  of  degree  n  -  1. 

We  also  define  a  forward  B-  spline  by 

(1.5)  Q„(X)  =  Mn(x-f)  «r‘  • 

Qn(x)  has  integer  knots,  is  positive  in  (0,  n)  and  zero  elsewhere. 

With  »jj  (t)  defined  by  (1.3)  we  define  • 

00 

(1.6)  $n(t)  =  %  »|<n(t  +  2itJ)  . 

J=-oo 

<j>  (t)  is  a  positive  cosine  polynomial  of  period  2u  and  order  [■— ]  -  1 

il  M 

that  can  be  explicitly  computed  from  the  expression 

0-7)  V*»  =  2  MnWe‘vt  =  En  Mn(v)eIvt 

V  =-00  \v\£  — 

By  Lemma  6  of  [9,  p.  180],  we  have 


(1. 8)  max  6  (t)  =  <J>  (0)  =  1,  min  <j>  (t)  =  4> (w)  >  0 

j  n  n  '  n  n 

By.  (1. 7),  we  find 


*2<t>  =  i 

♦3<‘>  =  1±fs± 


= 


2  +  cost 


*5<t)  = 


1 1 5  +  76  co s  t  +  cos  2t 


33  +  26  cos  t  +  cos  2t 

V*>  =  60 


#1183 


-11- 


We  shall  also  define  ^(t)  =  1  to  make  the  notation  in  what  follows 
more  convenient. 

In  the  Introduction,  we  referred  to  the  following 


Lemma  1.  If 


(1.9) 


then 


(1. 10) 


S(x)  c  Sn  n  L.(F) 
n  l 


2  I  s(vj i 


a.  id 


then 


(1.12) 


S(x)  c  Sn  n  Ljflf) 


Z  |S(»)|  <». 
o 


Proof.  Let  n=  2m.  We  reproduce  the  following  remark  of 
Louboutin  [6,  p.  l].  If  R^(x)  <  n^,  then  by  Markov's  inequal.  / 

(1.13)  max  |r*  (x)J  2  2k2  max  I R.  (x)  I  . 

fo.l]  *  [0,1]  k 

Let  now  P(x)  c  n_  ,  and  let 

CXOr  1 

X 

R2m(x)  =  /  P(t)dt  * 


Applying  (1. 13)  to  this  polynomial  of  degree  2m  we  find  that 

2  X  7  1 

max  |p(x)j  i  2(2mr  max  |  /  P{t)dt|  £  8rn  /  |P(x)|dx. 

[0,1]  [0,1]  0  0 

For  a  spline  function  S(x)  of  degree  2m  -  1  with  integer  knots,  we 


2- 


therefore  have 


«  rT* 

max  |$(x)|£8tn  /  |S{x)|dx. 

[v ,  v+1  ]  v 

Assuming  (1.9),  we  have 


00  00 

E  IS(v)l  2  £  max  I  S(x)  I  2  8m2||S(x)|L 

-00  .  00  [v,v+l] 

while  (1. 11)  implies 


oo  oo 

EiS(v)|  2  £  max  IS(x)|  2  8m2 1| S(x)|L  +  <  oo. 

0  0  [v,v+l]  > 

We  have  carried  through  the  proof  for  n  =  2mj  however,  we  ob¬ 
tain  the  same  result  for  n  =  2m  -  1  if  ws  replace  (1.13)  by 


max  |R'(x)|  S  2k2  max  Ir.(x)| 

[-1/2,  1/2]  *  [-1/2, 1/2]  * 


and  for  P(x)  t  n  define 
2m- 2 


R2m-l(x)  =  S  P<‘>dt  • 


This  completes  a  proof  of  Lemma  1. 


We  also*  need  to  know  Just  when  a  cardinal  spline  function  of 
degree  n-  1  is  in  Lj(F)  or  in  Lj(F+).  The  answers  are  given  by 

x  Suppose  S(x)  €  Sn  and  (1.4?  holds.  Then 


Lemma  2 


(1. 15) 


The  Inclusion 


Z  lc  I  <»  . 


v=-oo 


(1. 16) 


S(x)  «  L^JT) 


tolds  If  and  only  If 


(1. 17) 


Z  IcJ  <00  . 


Proof.  That  (1. 14)  is  equivalent  to  (1. 15)  is  a  special  case  of 
Theorem  12  of  [9,  p.  199],  and  is  hereby  established.  (1. 16)  follows 
from  (1. 17)  in  precisely  the  same  manner  as  (1. 14)  follows  from  (1. 15) 
in  the  proof  of  Theorem  12  of  [9,  p.  199].  We  now  start  with  the  as~ 
sumption  (1. 16)  and  wish  to  show  that  (1. 17)  holds.  This  is  derived 
from  the  previous  case  that  has  just  been  settled. Assuming  (1.4),  we 
consider  the  spline  function, 


(1. 18) 


Evidently 
(1. 19) 


S(x)  = 


rJSili .  i 

»=-[—j+ 1 


C  M  (x-v). 
v  n 


{S(x)  if  x  SO 

! 

0  if  x  =  -n  +  “, 


From  (1. 16)  we  conclude  that  S(x)  <  Lj(F)  and  the  first  part  of 
Lemma  2  shows  that  (1. 17)  holds. 


2.  Proof  of  Theorem  1  of  the  Introduction.  We  adapt  our  proof 
from  the  proofs  of  Theorem  2  of  [9]  and  Theorem  1  of  [13],  For  sim- 


-14- 


#1183 


plicity  we  write  M(x)  =  Mr(x)  and  =  H^t .  We  require  the  q.  f. 
(11)  to  be  exact  for  f(x)  =  M(x  -  j)  for  all  integers  j.  This  stipula¬ 
tion  gives  the  following  relations: 


(2.1)  /  M(x-j)e  <ix  =  ^n(t) e1 

-00 


w 

Z  H  M(v-J)  for  all  j 


v=-oo 


or,  since  M(v-j)  =  M(j-v)  as  can  be  seen  from  (1. 1),  we  have 


(2.2)  Z  M(J-v)Hy  =  4'n(t)e1 

v=-oo  * 


for  all  j 


To  Invert  this  convolution  transformation,  we  consider  the  positive 
cosine  polynomial  $  (t)  as  given  by  (1, 6}  and  the  expansion  of  Its 
reciprocal  in  a  Fourier  series: 


(2.  3) 


♦„<*> 


k=-co 


l  W*">  elkt 


Lemma  11  of  [9,  p.  187],  for  p  =  oo,  implies  that 
(2.4)  Hj  =  l  w£’  +n(t)ell,t 

is  a  bounded  linear  transformation  of  I  into  Itself,  whose  inverse 
is  given  by  (2.  2).  Since  the  sequence  {^(tje^*}  is  in  ,  we 

conclude  that  the  sequence  {H^  }  defined  by  (2,  4)  also  belongs  to 

t  ,  Since  wfn^  =  w^.  (9,  p.  182],  we  know  from  (2. 4)  that 
00  j- V  v-j 


(2.5,  H.  =  l  ,« ♦  I  .« 

V  =-00  i  V  =-00  1 


(n)  ei(v-j)t 


so  that  by  {2.  3)  we  get 
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<2. 6) 


H  -  M  .«» 

J  ‘  ♦„(*> 

The  sequence  {H. }  in  i  is  uniquely  defined,  by  (2. 4). 

J  00  ' 

A  proof  will  be  complete  as  soon  as  we  show  that  the  functional 


(2.?) 


w  .  w 

Rf  =  /  f(x)e  dx  -  Yi  H  £(y) 


v  =*<0 


with  the  coefficients  Hy  given  by  (2.6)  has  the  property  that 


(2. 8) 


Rf  =  0  if  f  c  S_  f|  L.(R)  . 

n  l 


Suppose  f(x)  is  such  a  function  and  let 


(2.9) 


f(x)  =  £  C  M(x-v) 


be  its  expansion  in  terms  of  the  central  B-  splines  of  degree  n  -  1. 
By  Lemma  2,  we  know  that  f(x) «  Lj(R)  implies  that 


(2. 10) 


Y  |C  I  <«o  . 


The  partial  sums 

(2.11) 


A 

yx)  =  Yj  CyM(x-v)  (r  «  0,1,2,...) 


v=-r 


have  the  additional  property  that 


(2.  12) 


yx)  =0  if  |x|  £  ~  +  r  . 


Since  f(x)  «  Sn  f)  Lj(R)  and  (2. 12)  holds,  we  conclude  that 

f  (x) «  S  n  L,(R).  Using  properties  of  the  functional  (2.7),  we  obtain 
r  n  l 


(2.13) 


f  yxje^dx  =  fj  HyW  . 


-16- 


#1183 


•v  v.-'  * v?* ;  si’-*'** ''  'r '  ">'  :  ¥‘r‘ 


p.  j^jfrrr5?rv^T^*<;  rM'T* 


Observing  that  each  ff(x)  is  dominated  by  the  function 

00  l 


£  Ic  |M(x-v) 


which  is  summable  on  R  by  Lemma  2  and  (2. 10),  we  see  that  on  let- 

j  .  ! 

ting  r  -*  oo,  the  relation  (2. 13)  becomes  the  desired  relation  1 

CO  .  00 

/  f(x)eixtdx  =  £  Hyf(v)  • 

-00  V=“ 00 

•  >  !  I 

This  completes  a  proof  of  Theorem  1. 

Substituting  the  coefficients  Hy  as  given  by  (2. 6)  gives  us  the, 

unique  q.  f.  of  Theorem  1  in  the  following  form: 

so  I*.  4*  (t)  oo  ,  . 

(2,14)  /  f(x)eiXtdx  =  -iL-  £  f(v)eivt+  Rf  . 

~?A  d)  (t)  V  =-00 

Yn  , 

,  i 

Suppose  now  that  f(x)  is  a  spline  function  of  degree  n  -  1  with 
knots  at  (v  +~)h,  for  all  integers  v ,  that  is  also  in  L^fR).  But 


then 


S(x)  =  f(xh) 


is  a  cardinal  sfiine  function  of  degree  n  -  1  for  the  step  1.  From 
Theorem  1  and  (2.  H)  we  have  1 

00  ,  .  ib  (t)  oo  ' 

/  f(xh)e  dx  =  ~rr  £  f(vh)e  v  .  , 

-co  ^n'  '  v=-oo 

,  i 

If  we  replace  in  the  integral  the  variable  x  by  x/h  arid  then  replace 

in  the  identity  t  by  th,  we  obtain 

00  ■  .  ib  (th)  oo  . 

(2.15)  /  f(x)elxldx  .  -“rj-  h  f  f(vh)eIwth  .  . 

-oo  mv  '  v  =-oo 
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y-y^yf  T>*TV^vt>-7 '•'7T'*  ?• 


■f*yry?r*i 


•r^t »~w -■■ istttvt 


w»vwr 

I 


If  f(x)  is  an  arbitrary  function,  then  this  is  no  longer  an  identity. 
However,  we  can  obtain  information  on  the  error  made  in  using  the 
approximation  (2. 15)  if  we  consider  n  even,  say  n  =  2m.  In  fact, 
in  §14  below  we  shall  prove  the  following 


Theorem  4.  Suppose  f(x)  c  C2m  n  L^R)  n  L2m(R)  and  ~ 

i  • 

is  a  natural  number.  Then  we  can  bound  I  Rf  |  as  given  in  the  q.  f. 

(2. 16)  /  f(x)eixtdx  =  T--7:uT  h  £  f(vh)elvth  +  Rf 
bZ 

(2.17)  |  Rf  I  £  4<~)2m  i|f(2m)L  for  all  rational  t*0 

I  IT  L>.y*>/ 

'  ln  i  Ijr  2 n. 

.ia  <-  h>  hK 

2  2  m 

In  the  theorem,  (R)  denotes  a  particular  choice  of  n  and 
p  for  the  set 

Ln(R)  =  {F(x)  :  F^n~^  absolutely  continuous,  F^eL(R)} 
P  i  P 

where  n  is  a  natural  number  and  1  £  p  =  «  .  The  set  Ln{R+)  is 

P 

defined  similarly. 


I 


I 
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-?;w.  '2*Z‘F%rf7Er^r*r&z+z^j& » V'”?* 


2L  THE  COSINE  AND  SINE  TRANSFORMS 


For  concreteness,  we  now  take  w(x)  =  cos  xt  and  consider  q.  f.  of 
the  form  (20).  We  shall  also  take  n  =  2m,  and  later  will  indicate  the 
modifications  necessary  for  different  derivative  data,  for  even  degree 
splines  (n  =  2m  -  1),  and  for  the  weight  function  sin  xt. 


3.  A  recurrence  relation.  For  simplicity,  we  write  ,  By 


B.'  .  \  We  want  to  construct  a  q.  f.  of  the  form 
h 1 


(3. 1)  /  f(x)cos  xt  dx  =  £  H  f(v)  +  £  B  f(2J~1)  (0)  +  Rf 

0  v=0  v  j=l  2J'1 


such  that 


(3.2)  |hJ  < K  for  all  Integers  v  ^0  and  some  K 

and  with  the  property  that 

(3.  3)  Rf  -  0  if  f  «  S,  n  L,  (R+). 

2tn  l 

We  do  this  by  enforcing  the  requirement  (3.  3)  for  an  appropriate  se¬ 
quence  of  elements  of  S,  D  L,  (H+).  The  sequence  we  require  is 

cm  i 

the  sequence  of  forward  B- splines  of  degree  2m  -  1 


(3.4) 


{Q(x  -  r)} 


(r  =  -2m+l,  -2m+2,  .,.) 


where,  by  substituting  n  =  2m  in  (1.  5)  we  have 


(3.  5) 


Q<*>  =  °2m(xi  =  (1^1)7  Eo  Ml'tfHx-i)* 


2  m- 1 


Since  Q(x)  has  support  in  (0,  2m), 
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VV^V'-'Va^ 


(3. 6)  Q(x  -  r)  =  0  outside  the  interval  (r,  r  +  2m), 

(r  =  -2rn  +  1,  -2m  +  2,  . . . ) 

Q(x)  «  n  on  any  interval  between  consecutive  knots  so  that 
2m-l 

Q(x)  j  S,  n  L.(JR+),  and  evidently  we  also  have 
m  m  x 

(3.7)  Q(x- r)  <  S2mn  L1(R+)  (r  = -2m+l,  -2m+2, . . . ). 
Substituting  f(x)  =  0(x  -  r)  in  (3. 1)  and  recalling  (3.  3)  and 

(3.  6),  we  have  the  sequence  of  relations 
r+2m 

(3.8)  /  0(X“r)cosxtdx  =  H  Q(-r)+H.Q(l-r)+. . . +H  0(2m+l) 

q  oi  r+zm-l 

+  E  ,Q(2J"1)(-r),  (r  =  -2m+l, -2m+2, . . .  ,-2,-1) 
j=l  J  1 

and 

r+2m 

(3. 9)  /  0(x-r)cosxtdx  =  Hr+1Q(1)  +  H.+2Q{2)+. . .  +Hr+2m_1Q<2m-l) 

(r  -  0, 1 ,  2, ... )  . 


4.  The  summation  of  certain  series.  We  shall  need  the  following 
lemma  which  deals  with  well-known  power  series. 


Lemma  3.  1*.  The  power  series 

oo 


(4.1) 


_  .  .  v»  k+1  v 

®k00  =  L  (v+1)  x 
v=0 


(k  =  0,1,2,...) 


has  the  sum 


(4.2) 


VX)  = 


pk(x) 


(1-x) 


k+2 
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where  Pk(x)  Is  a  monic  polynomial  of  degree  k,  with  lnteqer  coeffi¬ 
cients.  that  may  be  derived  from  the  recurrence  relation 

(4.  3)  Pk(x)  =  (l+kx)Pkl(x)  +  x(l-x)P'k  l(x),  with  PQ(x)  =  1  . 

2* .  The  power  series 


(4. 4) 


has  the  sum 


(4.5) 


*  (x)  =  2  (2v  +  l)k  xV 

v=0 


T  (x) 

=  k+l 
K  (1-x)  4 


where  Tk(x)  is  a  monic  polynomial  of  degree  k,  with  lnteqer  coeffi¬ 
cients.  that  may  be  derived  from  the  recurrence  relation 

(4.  6)  Tk(x)  =  [1  -  (2k-l)x]Tk  l(x)  +  2x(l-x)T*k  l(x)  with  TQ(x)  =  1. 

The  polynomials  Pk(x)  are  called  Euler- Frobenlus  polynomials  of  de¬ 
gree  k,  while  the  Tk(x)  are  called  mldpoli. :  Euler-Frobenius  poly¬ 
nomials  of  degree  k.  We  omit  the  easy  proof  by  induction  which  also 
furnishes  the  relations  (4.  3)  and  (4.  6).  We  find  readily  that 
Pj(x)  s  1+x 


(4.7) 


P2(x)  =  I  +  4x  +  x 

P3(x)  =  1  +  llx  +  llx2  +  x3 

P4(x)  =  1  +  26x  +  66x2  +  26x3  +  x4 
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(4.8) 


Tj{x)  =  1  +  x 
T_(x)  =  1  +  6x  +  x2 

C* 

T3(x)  =  1  +  23x  +  23x2  +  x3 

T,tx)  =  1  +  76x  +  230x2  +  76x3  +  x4 
4 


and  so  on. 


The  form  of  the  relations  (3. 8)  and  (3. 9)  suggest  the  use  of  gen¬ 
erating  functions  for  the  determination  of  the  and  the  .  The 

righthand  side  of  (3.  8)  and  (3. 9)  is  equal  to  the  coefficient  of 

r+2m- 1  . 
x  in 

oo  2m- 2  2m—2  m-1 

(4.9)  (£  Hxv)(  j  0(2m-l-v)xV)  +  E  [  Eb 2\  1  (2m-l-v)]x 

v=0  V  v=0  v=0  J=1 

In  order  to  simplify  the  two  polynomials  in  (4. 9),  we  note  that 


(4. 10)  0(k)(x)  =  (-l)kQ(k)(2m-x) 


(k  *0,1,.,.,  2m- 2) 


as  can  be  verified  from  (3.  5).  With  this  substitution  and  the  inter¬ 
change  of  the  order  of  summation  in  the  second  polynomial,  (4.9)  be¬ 


comes 


2m- 2 


2  m- 2 


(4.11)  (E  HxV)(  E  Q(v+l)xV)-E  B  (  E  Q  2,-1  (v+l)x  ) 


We  need  the  following  result  that  is  perhaps  of  independent 


Interest; 


Theorem  5.  The  following  identities  hold: 
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li. 


2m~2  .  .  (l-x)'P,  (x) 

(4.12)1*.  J  QU,(v+l)xV  =  "  0=0,1,...,  2m-2j 

m=l,  2, . . . ) 


m- 1 


.{2  m- 1). 


,2m- 1 


(4.13)2*.  £  Q  *'(v+l-0)x’  =  (1-x)—  *  (m  =  1, 2, . . . ) 


v=0 

(111I3-  2l?2iii,))  |L(1  nk”  ,1.2m- 2-1  (1-x)‘T2m-2-i,x) 

v=0  2m‘1  *?  (2m-2-J)I 

(J=0, 2m- 3;  m=l,2, ...) 


2  m- 2 


4*.  2  M 


v=0 


(2m-2) 
2m- 1 


(v+l-m-0)x  =  (1-x) 


2  m- 2 


( m= 1 ,  2, . . . )  . 


We  note  that  because  of  (1.  5)  we  could  also  write  (4. 12)  in  the  form 

2,£2  m  v  <>-x)JP  j  i(x> 

(4.15)  T.  M"'(v+i-rn)x  =  - SatA-J - 

w  2m 


v=0 


(2m-l-j)! 


0=0, 1,. . .,  2m-2j  m=l,  2, ...) 


Proof.  First  we  show  (4. 12).  Let 
.  2m-  2 

(4. 16)  (l-x)JP  (x)  =  J  A  xV  • 

2m- 2- j  v 

From  (4. 1)  and  (4.  2),  for  k  =  2m  -  2  -  j,  we  find  that 


00 


(1’x)  P2m-  2-J(x) 


=  (l-x)2m£  (Ml)2ra-1->xk 
k=0 

2m  ...  oo 

=  >:  (-i)‘(2Im)X1- 1  <k+i)2 

1=0  k=-oe 


2m-l-j  k 
x 


so  that 
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.  V/  ..1,2m...  , ,.2m-l-J  v  /  ..1,2m.,  ,,,.2m-l-J 
(4.17)  A  =  2j  (-1)  (  .  )(k+l)  =  2j  (-.\)  (  )(v-i+l) 

V  i+k=v  +  i=0  1  + 

On  the  other  hand,  by  differentiating  (3.  5)  j  times,  we  obtain 


t4-,8»  q(i>m  -  &  <-i'i(2r^i.r1-1 . 


and  then 


(4.19)  (2m-l-J)!  Q  (v+1)  =  £  C-l)  (  !")(v-i+ir 

i=0  r 


which  is  identical  to  (4. 17)  so  that  (4. 12)  follows. 

To  prove  (4. 13),  we  substitute  j  =  2m  -  1  in  (4. 18)  and  get 

(4.20)  Q(2m_1)(x)  =  J  (-uV^'Mx-  l)°  . 

1=0  + 

p(2m-l)^j  lg  a  step  function  so  that  upon  substitution  of  x  =  (v+1-0) 
for  v  =  0, 1, . . . ,  2m-l,  (4. 20)  becomes 

(4.21)  Q(2ni“I) (v+1-0)  =  ^HjVjV+l-i-O)?  *  £(- l/(2m) 

J  <4  •  T  J  A  • 


(v  =  0, 1 , . . . ,  2  m- 1 )  , 


But  an  easy  induction  shows 


(4.22)  £  (-UVI")  =  (-1)V(^“1)  (v=  0,  l,...,2m-l) 

i=0 

so  that  by  (4.  21)  and  (4.  22),  we  have 


2m-l  , .  2m-l  ,  _  , 

l  Q(2m-1Vu-0>*1'=  l  (-1)V  (  V  =  (l-x)" 


The  even  degree  spline  case  3° ,  4*  is  proved  in  the  same  manner 
as  (4.  12),  (4.  13)  so  we  can  omit  this  proof  and  the  theorem  follows. 
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With  the  substitution  of  (4. 12),  (4. 11)  becomes 


»  P-  ,(x)  m-1 

(4.23,  (Zh/I-^JT  *2 


,"y 1  „  (x~]>  1  f2m-2l-l(1Ii 

2M  (2m-2J)f 


We  now  turn  our  attention  to  the  left  side  of  relations  (3. 8)  and 
(3. 9)  and  define 


r+2m 


Fr+2m-l  =  f  O(x-r)cos  xt  dx  (r  =-2m+l,-2m+2, . . .  ,-l,  0) 


(4.  24) 


r+2m 


r+2m-l 


*  /  Q(x-r)cos  xt  dx  (r=l,2, . ..). 


If  we  integrate  the  righthand  side  of  (4.  24)  for  (r  =  -2m+l,  ...,-1,0) 
by  parts  2m- 1  times,  we  obtain 


(4.  25)  Fr+2m- 1  =  f'tC^x”-r)s*n  xt  +  •“  Q'  (x-r)cos  xt  -  ~  Qn(x-r)sin  xt 

,  .  ..m-1  ..  .  r+2m 

-  *7  Q"'(x-r)cos  xt  +  . . .  +  i^— pO*  Vr)sin  xt] 

t4  tZ  1  0 


..m-1  r+2m 

—  f  0 

t2*-1  0 


(2m- 1) 


(x-r)sin  xt  dx  . 


Since  Q^2m~1\x)  is  a  step  function,  we  break  up  the  interval  of  inte¬ 
gration  and  writ''  the  integral  in  (4.  25)  as 

(4.26)  /  Q(2tn“1)(X_r)Sin  xt  dx  =  f  1Q(2m> l)(x-r)sin  xt  dx 

0  1=0  1 


r+2m-l  ,,  ..  f+1 

£  q (2m-l )(|+i_r-o)/  sin  xt  dx 

1=0  I 


After  we  integrate  and  collect  terms,  we  obtain  (4.  26)  in  the  form 
#1183  -25- 


r+2m-l 


1  *  1  Mill  *  |  4 

(4. 27)  -  £  Q*  '(#H-r-0)  [cos(l+l)t  -  cos  ft] 

1=0 


r+2ra-l 


=  -“{-O(2m-1>(l-r-0)  +  j  [Q r~0)~0  <2m“1>{jC+l-r-  0)J 


1=1 


•  cos  it  +  Q^2m”^(2m~0)cos(r+2m)t}  . 


By  (4. 21)  we  find 


(4.28)  Q(2m“1)(|-r-0)  -  Q*Zmrl){M-r-Q)  =  -(-l)1"^^) 


and 
(4.  29) 


Q(2m“1){2ir*-0)  =  £  (-l)1^  ) 

1=0 


=  -  l 


so  that  upon  substitution  in  (4. 27),  we  have 
r+2m 

f 

0 


(4.30)  /  Q^^x-rjsinxtdx 


1  n(2m-l).  .  ,  1  r+v>m,  . .f-r , 2m 

=  ~Q'  '(l-r-0)+-  h  (-1)  ( jj_|-)cos  It  • 

If  we  let  i  =  2m  +  r  -  I  in  the  sum  in  (4.  30)  and  substitute  the  result 
in  (4. 25),  (4.  25)  becomes 


(4. 31)  F, 


r+2m-l 


,m 


\  Q’(-r)  Q'"(-r)+. . .  Q(2m'3)(-r) 


m  r+2m-l 


+14ir  W»>  i  (ii)1(2r)cos(r+2m-i,t 

tZm  t2m  uo  1 


(r  =  -2m+l,  ...,-1,0). 
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Let  J  =  r+2m-l  and  multiply  each  side  of  (4. 31)  by  (-l)mt2m  to  get 


(4.32)  (-l)mt2raFj  = -[(-l)mt2m“2Q,(2m-l-j)-(-l)mt2m"40,"(2m-l-JH 

...+  t2Q(2m“3)(2m-l-j)]  +  Q(2m“1)(2m-j-0) 

+  ^  (-dVJVos  <j+l-i)t  (1*0,1,..,.  2m~l) . 

1=0 

In  view  of  the  relations  (4. 10),  we  may  write 
(4.  33)  (-l)mt2mF  =  £  (-l)i(2m)cos(j+l-i)t  - 

1  i=o  1 

-Q(2m-l)(j+i-0)+t2Q(2m-3)^+1)_  ^  .+(-i)m-lt2m-4QIU(j+1j 

+  (-l)mt2m"2Q'(j+l)  (j  =  0,  l,...,2m-l). 


If  we  consider  (4. 24)  for  r  =  l,  2, . . .  and  again  integrate  by  parts 
2m- 1  times,  we  obtain  (4.  25)  with  the  lower  limit  of  integration  r 
instead  of  0.  Upon  evaluation,  the  square  bracket  in  (4.  25)  is  0  and 
we  would  get,  by  using  (4.  25),  (4.  26)  and  (4.  27)  that 


,  ,.m  f  .  r+2m-l 

(4  34)  F  -  inil —  J  JL  V  0(2m-l) 

(  *  }  r+2m-l  ~  2m-l  )  t  Q 

t  [  I=r 


(i+l-r-0)  [cos(I+l)t-cosj 


(r  =  1,2,...). 

Following  the  same  steps  we  did  before,  and  noting  by  (4,  21)  that 


allows  us  to  use  (4.  30)  to  write 


,  ,.ro  r+2m  ,  ,, 

(4-J5>  l  <1>S  *  • 

Letting  1  =  2m  +  r  -  i  in  the  sum  In  (4.  35)  and  then  J  =  r  +  2m  -  1, 
we  obtain  upon  multiplication  by  (-l)mt2m 


(4.36)  (-l)mt2mF.  =  2  (-l)1(Z|n)cos(j+l-i)t  (J=2m,  2m+l, . . . ) 

1  1=0  1 


,1.2m 


Let  us  sum  the  series 

(4.37)  (-l)mt2m£  FxJ 

j=0  1 

from  (4.  33)  and  (4.  36)  we  find  that 


(4.38)  (-l)mt2m  E  F4xi  =  (  Z  t  Tj  (- 1  )*( 2!”)cos(j+l-i)t]x* 
J=0  i  [  j=0  1=0 


Q(2m-!)(J+i  0)XJ 


»  2m  .  ,  A  2m-l 

+  E  fZ  (“frV. "Wo+l-DtjA-  E 

J=2m  1=0  )=0 


2  m- 2  .  _  2  m- 2 

+  t2  E  Q  2  3  o+i)x“. .  *  +<- i)mt2m~2  E  Q'o+Dx . 


To  simplify  the  term  in  curly  brackets  on  the  right  side  of  (4.  38),  we 


define 


(4.  39) 


t(x)  =  E  [cos{v+l)t]x* 
v=0 


and  note  that 

7  2m  .  .  °° 

(4.40)  (1-x)  m-r(x)  =  E  (”U  (  I”)*  *  E  [cos(v+l)t]xV 

1=0  v=0 

co  2m 

=  E  l  Z  M)  (  ["icosCv+Dtjx1 

j=0  i+v=j 

1=0.  v  SO 


29 


^  ^  00 
=  Z  l  ZmjV 1m)cos(j+i-l)t]xJ  +  2  [  2  (-l^^f'jcosO+l-Dtjx5. 

J=0  1=0  j=2m  1=0 

Using  (4.40)  and  (4. 12)  and  (4. 13)  of  Theorem  5,  we  may  therefore 
write  (4.  38)  as 


(4.41)  (-l)mt2m  £  FjX*  =  (l-x)2mr(x)  -  [(l-x)2™"1 


J=0 


~j  PjixJd-x)2"1”3  P3(x)(l-x)2m~5 


+  (-1) 


m-1  t 


2  m~  2 


(2m- 2)!  2  m- 3 


P,_  Jx)(l-x))  . 


00  1 

Equating  the  relation  (4.  23)  and  F  x1  as  determined  from  (4.  41), 


J=0 


we  see  that  we  require 

.m 


(4.42)  {(x-l)2mT(x)+(x-l)2m'1-|rP1(x)(x-l)2m'3  + 


i  .2  m- 2 

+  *~1)  (2m- 2)!  P2m-3W<X_1^ 


oo  P,  ,(x)  m-1 

.  y  B 


-<l  HxV 


v=0  V  '  ( 2  m—  1 ) !  2J-1  (2m-2j) ! 


5.  Determination  of  the  coefficients  H  =H^2m^.  B.  =  b!2!^ 

■»  »  V  Vy  t  j  j,t 

00 

Solving  (4.  42)  for  £}  H  xv  gives  the  final  relation 

v=0  V 

(5.1)  f  H/  =  p2,t?'j'(x)  fcitr  t(x-l)2tnT  (x)4-(x-l)2m~1  + 
v=0  2m-2'  ;  t 
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R£l  i  t2j 

+  1  (-i)J  TrrrrP 


i=i 


(2j)l  2j-l 


(x)(x-l) 


2m-2j-l 


”^1  (x-l)  J  p2tn-21-? 


1  -  S  Vi 

j=l  J 


(2m-2J)f 


Our  derivation  of  (5.1)  establishes  the  following 


Proposition  1.  The  coefficients  H  =  H^2m^.  B_ .  .  =  Bi2m?  . 
-___c -  -  v  v,  t  *  2j- 1  2j-l,t 


of  the  most  general  functional 

oo  eo  IDrl 

(5.2)  Rf  =  /f(x)cos  xt  dx  -  Yt  H  f(v)  -  j)  B  f  2^*X*(0) 

0  v=0  V  j=l 

vanishing  for  the  functions 

(5.  3)  Q2m*X  “  r)  (r  =  -2m+l,  -2m+2,  . . . ) 

are  the  expansion  coefficients  of  the  rational  function  (5. 1)  where  the 

B,  (j  =  1 . m-1)  are  chosen  arbitrarily. 

ZJ-1 

We  want  to  investigate  the  functionals  (5. 2)  further,  and,  in  par¬ 
ticular,  determine  the  unique  functional  having  bounded  coefficients. 

Let  R,  (x)  denote  the  right  side  of  (5. 1),  where  the  B_.  ,  (j=l, .. . ,  m- 
2m  2j-i 

are  as  yet  undetermined.  To  use  R,  (x)  effectively,  we  need  informa- 

2m 

tion  on  its  poles.  To  this  end,  in  view  of  (4. 12)  of  Theorem  5  for  j  =  0, 

we  may  write  P,  ,(x)  in  terms  of  the  central  B-spline  (1.  5)  as 
2  m- 2 

2m-2 

(5.4)  P2m-2*X)  =  f2™-1)*  Z  M2m(v-m+l)xV. 

v=0 

By  Lemma  8  of  [9,  P.  182]  we  know  that  this  reciprocal  polynomial  has 
only  simple  and  negative  2eros  so  that  we  may  label  them  to  satisfy 


the  conditions 


(5.  5) 
and 
(5. 6) 


X2m-2  <#"  <  \n  <~1  *  Xm-1  <•••  <  xj  <0 


XlX2m-2  “  X2X2m-3 


•  •  •  ~  i  ^  ~  1 

m- 1  m 


From  the  form  of  t(x)  as  given  by  (4.  39),  we  note  that  r(x) 

i  i  1 

converges  for  |x|  <1.  Observing  that  for  R_  (x)  the  poles 

CX!i 

X. .....  X  ,  are  inside  the  unit  circle  while  X  .....  X  _  are  out- 

l  m-i  m  m -c 

.  •  i 

side,  in  view  of  (5.  5),  we  shall  have  the  coefficients  bounded, 

if  and  only  if  the  coefficients  B  (j  =  1, . . . ,  m-1)  can  be  chosen 

so  that  the  m-1  poles  X-.....X  ,  of  R_  (x)  have  vanishing  re si- 

l  m-i  cm  ■ 

dues.  By  (5. 1)  this  will  occur  if  and  only  if  the  B.'  satisfy  the 


equations 
m-1 


(X-l)2,“lP,m  ..  .  (X  ) 
v _  2m-zi-l  v 


(5,7)  (2m- 2  j)! 


m 


_  til 

^2m 


[(X  -l)2mr(X  )+(\-l)2m’1+  2H)J  — 
v  v  v  )=1 


(2J)! 


‘  P2j-l(XvKV1,2m2J"lj  ‘  (v  =  l,...,m-l)  . 

In  the  system  (5.7),  we  have  m-1  equations  in  the  ni-1  unknowns 
B_.  (j=l, . . . ,  m-1).  To  show  that  the  system  is  nonsingular,  it  is 

sufficient  to  show  that  the  determinant 


(5.8) 


(\-l)2j  *P 

U  I  =  - 

*  vj*  1  (2m-2j)f 


2.rni2h).(^).  I  /  0 

1  \  t  *  f  U 


(v  =  1,2,...,  m~  1)  j=l,..»,  m—  1 
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In  order  to  accomplish  this,  we  shall  consider  a  related  problem, 
a  special  case  of  which  will  solve  our  problem.  Let 

(5.9)  =  {S(x)  :  S(x)  c  S2m,  S(v)  =  0  for  all  integers  v}  . 

In  (9,  p.  194]  Schoenberg  shows  that  every  element  of  s!?  admits  a 

cm 

unique  representation 

2  m-  2 

(5.10)  ,  S(x)  =  £  aiA(x> 

k=l  K  K 

for  appropriate  values  of  the  coefficients  a. ,  where  the  S  (x)  are 

K  K 

the  so-called  elgenspllnes  of  the  class  S®  and  are  defined  in  terms 

Zm 

of  the  zeros  (5.  5)  by 

00 

(5.11)  Sfc(x)  =  S2m(x;  y  r  S  xjcM2m(x",)  (k=1>  2, . . . ,  2m-2)  . 

J=-00 

In  [9,  §9]  Schoenberg  proved  a  Theorem  11,  a  special  case  of  which 

! 

asserts  the  following: 

(5.  12)  S(x)  <  S2m  fl  L®(F)  for  some  s  =  0,  1, . . . ,  2m- 1 

implies  that 

(5.  13)  S(x)  =  0  for  all  x. 

The  first  half  of  the  proof  actually  establishes  the  following: 

Every.  S(x)  «  S2m  n  L  ®(F+)  for  some  s  =  0,  1, . . . ,  2m- 1  may 
be  uniquely  represented  in  the  form 

m- 1 

(5.14)-  S(x)  =  £  a.S.(x) 

k=l  K  K 

for  appropriate  values  of  the  coefficients  . 


-32- 


#1183 


In  particular,  the  S^{x)  for  k  -  1,  2, . . . ,  m- 1  are  linearly  inde¬ 
pendent. 

We  determine  the  set 
(5.  15)  I  C  {1, 2, . . . ,  m- 1 } 

the  null  set  also  being  allowed,  its  complement 

1°  =  {l,  2, . . . ,  m- 1 }  -  I 

c  * 

and  in  terms  of  I  ,  the  set 

(5.  16)  I*  =  {2m-2-i  :  i  c  IC>  . 

Notice  that  while  lisa  subset  of  {l,  2, .  . . ,  m- 1 },  the  new  set  I* 
is  a  subset  of  {m,  m+1, . . . ,  2m+2  }  so  that  the  intersection  I  n  I*  is 
empty.  The  particular  set  we  will  be  concerned  with  in  connection 
with  the  system  (5.  7)  is  I  U  I'=  {l,  3,  5, ... ,  2m-3  }  . 

Suppose  that  S(x)  is  of  the  form 

m-1 

(5.17)  S(x)  =  Yj  <Wx). 

k=l  K  K 

We  want  to  be  able  to  choose  the  in  such  a  way  that 

(5.  18)  S(i,(0)  =  yQ(i)  i  €  I  U  I' 

where  the  righthand  side  has  arbitrarily  prescribed  values.  In  other 

words,  we  want  the  a^  to  satisfy 

(5.19)  £  aiX1)(0>  =  yi1)  iclUP. 

k=l  *  K  0 

This  then  is  the  related  problem,  a  solution  to  which  will  enable  us  to 
show  (5.  8). 
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In  order  to  show  the  existence  and  uniqueness  of  the  a^,  we  con¬ 
sider  the  corresponding  homogeneous  system  and  prove  the  following 

0  &  + 

Lemma  4.  If  S(x)  c  H  L^  (R  )  for  some  s  =  0, 1, . . . ,  2m- 1 

and 

(5.20)  S^(0)  =  0  if  1  s  I  U  I* 

then 

(5.21)  S(x)  =0  for  all  x. 

Remark.  This  is  a  special  case  of  a  result  of  Schoenberg  [10, 
Lemma  2]  concerning  a  finite  interval  that  we  have  extended  to  the 
infinite  interval  (0, «).  We  follow  his  proof  which  was  in  turn  Dased 
originally  on  a  proof  of  Greville  [4,  p.  4], 

Proof  of  Lemma  4.  Suppose  (5.  14)  is  the  canonical  representa- 
0  s  + 

tion  of  S(x)  «  S^^  n  L j  (F  )  for  some  s  =  0,  1, . . . ,  m- 1.  Note  that, 
since 

(5.22)  sj^  c  Lj(R+)  for  (s  =  0,  1, . . . ,  2m- 1)  and  k  =  1,  2, . . . ,  m- 1 
by  the  nature  of  the  representation  (5.  14),  we  also  have 

(5.23)  S(s)(x)  «  LJ(R+)  for  s  =  0,  1, . . . ,  2m- 1 . 

We  let 

(5.24)  12  =  /  [S(m)(x)]2dx  =  lim  f  [S(m)(x)]2dx 

0  b  **oo  0 

and  wish  to  show  that 

{5.  25)  n  =  0. 
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r  fm\  / m \ 

We  write  n  =  /  S'  'S'  ;dx  and  integrate  by  parts  successively  to 
obtain 

(5.26)  fi  =  lim  /  S(m,dS(m"1)  =  -  lim  /  S(m+1  Wm"2>  =  . . . 

0  0 

=  +lim  /  S(P)dS(tt) 

0 

where  the  integers  a  and  p  satisfy  the  conditions 
(5.  27)  0  £  <?  i  m-1,  mips  2m- 1,  a  +  p  =  2m-l . 


Notice  that  we  have  written  (5.  26)  as  if  all  the  ’’finite  parts"  drop  out 
at  each  end  of  the  interval  of  integration  and  at  each  step  of  the  suc¬ 
cessive  integrations.  That  this  is  indeed  the  case  follows  thus:  for 
each  pair  of  numbers  (a,  p)  satisfying  (5.  27)  we  have 

(5.28)  S(a)(0)S^(0)  =  0 

because  either  a  «  I  or  p  «  I’  so  that  (5.  20)  implies  that  (5.  28) 
holds.  For  each  ( a ,  p)  we  also  have 

lim  S(or)(b)S((3)(b)  =  0 

b  -*  oo 

in  view  of  (5.  23). 

The  integrations  can  be  continued  all  the  way  down  until  we 
reach 

b  .  [bl  Xv+1  , 

(5.  29)  n  =  lim  /  S(2m"  JdS  =  lim  /  S(2m" ^(xJS'Mdx 

0  b-*  co  v=0  X 

v 

where  we  have  written  x^  =  v  for  v  =  0, 1, . . . ,  [b]  and  x^j+j  =  b  • 

Since  S^2m  ^(x)  is  a  step  function,  the  integrals  in  the  sum  (5.29) 

vanish  if  0  =  v  <  [b],  since  S(x)  e  .  There  remains  to  show  that  also 

’  2m 
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(5.30)  lim  J  S^Zm  ^(x)S'(x)dx  =  0. 
b-*-»  [b] 

If  b  is  an  integer,  (5. 30)  is  true.  If  b  is  not  an  integer,  since 
S<2m~l\x)  is  a  step  function, 

(5.31)  /  s^Zm”1^(x)S,(x)dx  =  S^Zm**1^(b)  /  S'(x)dx 

[b]  [b] 

=  S{2m"1,(bHS(b)  -  S([b])]. 

S(x)  <  S®  implies  S([bl)  =  0  and,  upon  letting  b-*  <x>  on  the  right 
zm 

side  of  (5.  31),  we  obtain  (5.  30)  by  virtue  of  (5. 23). 

We  have  just  established  (5.  25),  and  therefore  that 

(5.32)  S(x)  «  . 

But  S(x)  c  S®  implies  that  S(x)  s  0.  This  completes  a  proof  of 
2m 

Lemma  4. 


In  view  of  Lemma  4,  S(x)  as  given  by  (5. 17)  for  the  homogeneous 


system 
(5.  33) 


m-1  ... 

I  «ks  'l(0]  =  0 

k=l 


i  c  IUI' 


must  vanish  for  all  x.  Then,  since  the  S^x)  for  k  =  1, ...,m-l  are 
linearly  independent,  we  must  have  a^  =  0,  k  -  1, . . . ,  m-1  .  This 
shows  the  existence  and  uniqueness  of  the  a^  for  the  system  (5. 19), 
so  that  we  must  have  that  the  determinant 


(5.34)  |Aikl  =  |sj.1,(0)|#0  (it  I  U  I*,  k  =1,2,...,  m-1). 
We  claim  that 
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(5.  35) 


S^W) 


<V‘>  P?.,W.-.<V 

(2m-l-i)f 


m- 1 


(1=0, 1, . . . ,  2m-2;  k=l, . . . ,  m-1). 

Indeed,  if  we  differentiate  (5. 11)  i  times  and  substitute  x  =  0,  we 

* 

obtain 


=  l 


(5.  36) 

n  . 

)=-00 

M„  (x)  has  support  in  (-m,  m)  and 
cm 


(5.  37) 


as  can  be  seen  from  (4. 10),  so  that  we  can  write  (5.  36)  in  the  form 

(5.28)  S(1)(0)  =  (-1)1  t  4.  M2m(J)  * 

j=-(m~l) 

If  we  let  v  =  j  +  m  -  1  in  the  sum  in  (5.  38),  we  obtain 

«5.3„ 

\  v*° 

Then  (4. 15)  of  Theorem  5  establishes  our  claim  (5.  35). 

th  rn  *  1 

If  we  multiply  the  k  column  of  the  determinant  in  {5.  34)  by  X^ 

k  =  1, . . . ,  m-1  and  take  the  transpose  of  this  resulting  determinant, 

we  have  by  (5.  34)  that  the  determinant 

(X.-l)‘p,m 

(S.  40)  I-  (Zmlw)!  "  ^le,  l)  /  0  (k=l,2, ....  m-1,  lcIU  I'). 

If  we  consider  the  special  case  I  U  I'  =  {1,  3,  5, ... ,  2m- 3},  then  the 
determinant 
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,5  41)  |  (Vli‘‘W2-l(V| 

(5*41)  1  (2m-l-i)! 


(k,l) 


/  0 


(k=l,  2, . . . ,  m-lj  i=l,  3,  5, ... ,  2m-3). 
This  is  precisely  the  relation  (5. 8),  so  that  we  have  established  the 
first  part  of  the  following 


Theorem  6,  1*.  Among  all  functionals 


00 


00 


(5.  2)  Rf  =  /  f(x)cos  xt  dx  **  H^}f(v)  -  £  B™ 

0  v=0  *  j=l  J  * 

vanishing  for  the  sequence  of  spline  functions 


m-1 


(2m)  f(2J-l)(0) 


(5.  3) 


(r  =  -2m+l,  -2m+2,  . . . ) 


there  is  a  unique  one  such  that  the  sequence  (H^  }  is  bounded. 


2*.  In  factr  this  unique  functional  Rf  can  be  given  explicitly  by 


00 


00 


+,m(t)  f 

(5.42)  Rf  =  J  f(x)cosxtdx-  <  ~f(0)  +  f(v)cos  vt 

0  *2m(t)  I  2  v=l 


,  V  (-l)Jf  4>2j(t)4j2m-2J(t)  ]  f(2J-l) 
U  ill1”  j.  /t\  J1 
J=1  t 


&  1‘  *2m(t) 


(0) 


(where  d>  (t)  and  <j>  (t)  are  defined  in  (1.  3)  and  (1.  6),  respectively) 
n  n 

for  m  =  1,  2,  3,  4. 


•Proof  of  2*  of  Theorem  6.  We  observe  that  the  functional  Rf  of 
(5.  42)  is  of  the  proper  form  (5.  2)  where 
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(5.  43) 


=  If  2|  •  •  •  )  9 


» 


„  »  ^ 
0  =  2 


H 


By  (1. 3)  and  (1. 8), 


cos  vt 


(t) 

(t) 


(5.44) 


!h  i  * 


sin  t/2 


t/2 


2m _ 1 


*2m(’') 


*ZmM 


<  00 


(v  =  0, 1,  2, ... ) 


so  the  sequence  {H^}  as  given  by  (5.43)  is  bounded.  Once  we 
show  that  Rf  as  given  by  (5.42)  vanishes  for  the  sequence  of  splines 
(5.  3),  the  unicity  established  in  part  1°  will  establish  part  2° . 

We  accomplish  this  in  two  steps,  one  for  (r  =  0, 1, 2, . . . )  and 
the  other  for  (r  -  -2m+l,  -2m+2,  . . . ,  -1).  We  remark  that  we  prove 
the  first  case  for  general  m  but  the  latter  only  for  the  special  cases 
of  m  =  1,  2,  3,4.  So  far,  the  latter  general  case  still  eludes  us;  it  is 
a  matter  of  showing  the  validity  of  one  necessary  identity.  This  same 
dilemma  prevents  us  from  claiming  explicit  versions  of  other  q.  f,  to 
follow  as  well  as  (5.42)  for  general  m. 

We  first  show  that  Rf  =  0  for  f(x)  =  Q,  (x-r)  for  r  =  0, 1,  2, . . .  . 
Since  (1.2)  holds,  we  get 

/  Mn(x-j)eixtdx  =  . 

-00 

Taking  real  and  imaginary  parts,  we  find 

00 

(5.45)  f  M  (x-,j)cos  xt  =  4/  (t)  cos  Jt 

n  n 

-00 

and 
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60 

(5. 46$  /  Mn(x-j)  sin  Jt  =  ^(t)  sin  jt . 


-00 


By  (1.  5)  and  (5.45)  for  n  =  2m  and  since  M,  (x)  =  0  if  |x|>  m, 

2m 

we  obtain 


=  *2m(t,cos  (m+r)t  (r  =  0, 1, 2, . . . ). 

Since  f(2J  ^(0)  =  0  for  f(x)  =  Q2m(x-r)  and  r  =  0, 1,  2, . . .  we  need 
only  show  that 


r+2m-l 

(5.48)  4>  (t)cos(m+r)t  =  Y.  M  (w-m-r)cos  »t 

Zm  v=m  2m 


(r  =  0»  lj  2^  . . . ) , 


By  (1.  7)  we  have 


m-1 


=  JL>  M2™<k,cos  kt 


which  upon  multiplying  by  cos(m+r)t  and  letting  v  -  m+r+k  becomes 


r+2m-l 


(5,49)  4>2m(t)cos(m+r)t  =  ^  M?m(v-rn-r)cos(v-m-r)t*  cos(m+r)t. 


v=r+l 


2  m' 


By  using  the  identity 

oos  at  cos  bt  =  cos(a+b)t  +  sin  at  sin  bt 

we  get 

(5.  50)  cos(v-m-r)t  cos(m+r)t  =  cos  vt  +  sin(v-m-r)t  sin(m+r)t. 
But 
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r+2m-l 

T.  M,  (v-m-r)sin(v-m-r)t- sin(m+r)t  = 

. »  2m 
v=r+l 

m- 1 

=  sin(m+r)t  Y.  M.  {k)sin  kt  =  0 
k=-(m-l)  2m 

since  M2m{-k)  =  ^2m^’  S0  tkat  upon  su^stitutin9  (5.  50)  into  (5.  49) 
we  obtain  (5.48)  as  we  wished. 

What  remains  then  is  the  case  r  =  -2m+l,  -2m+2, . . .  ,-l.  By 
(1.  5)  we  can  write  (5.3)  as  M,  (x-m-r)  and  need  to  show  that 
Rf  =  0  for 


or  for 


(5.51) 


M2m(x-m-r) 


(r  =  -2m+l,  -2m+2,  . . .,  -1) 


(j  =  -m+l,  -m+2,  . ..,  m-1). 


By  the  symmetry  of  M2m(x)  and  (5.45)  for  n  =  2m,  we  find 

00  oo 

(5.  52)  f  M2m(x-j)cos  xt  dx  =  ^2m(t)  cos  jt-  /  M2m(x+j)cos  xt  d> 


and 


00 

(5.53)  /  M  (x)  cos  xt  dx  =  “  (t). 

-  2m  2  2m 


2  2m 
u 

Since  M(Zk"1)(-j)  = -M(2k~  U(J)  and  M(2k“1)(0)  *  0  for  (k=l,2,..., 
m-1),  and  by  the  previous  case  for  (r  =  0, 1, 2, . . . ),  we  need  only 
show  that  Rf  r  o  for  (5.  51)  for  (j  =  -m+l,-m+2, . . ,  ,-l)  .  For 
m  =  1,  2,  3,  4  this  is  just  a  matter  of  computation.  For  instance,  for 
m  =  2,  the  cubic  case,  we  need  show  only  that 

r  *4(t)  1  l  <Mt)4,2(t) 

(5.54)  /M4(x+1)cos  xtdx  ^M4(l)}-”[1- - 


♦4(« 
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By  integrations  by  parts,  the  left  side  becomes 

1  t2 

-[cost-  1+  — 1 

and  the  right  side,  upon  substituting  M^d)  =  1/6,  =  -1/2, 

2  co  s  t 

4»2(t)  =  1,  <j>4(t)  = - - -  and  using  (1.3)  and  trigonometric 

manipulations,  agrees.  This  establishes  Theorem  6. 


We  can  now  prove  the  following 


Theorem  7.  Among  all  q.  f.  of  the  form 


CO 


m-1 


/  f(*)cos  xt  dx  =  l  H'2”V)  ♦  l  B<2“  +  Rf 

0  v=0  *  J=0  *  ' 

where  the  H^2t^  satisfy  the  condition 

IH^^I  <  K  for  fixed  t,  for  all  integer  v  =  0  and  some  K, 

there  is  a  unique  q.  f.  with  the  property  of  being  exact,  i.  e. ,  Rf  =  0, 

4. 

whenever  f(x)  t  Lj(R  ).  This  unique  q.  f.  Is  given  by  (5.  42) 

for  m  =  1,  2,  3,  4, 


Proof.  The  proof  is  modeled  after  the  proof  of  Theorem  1.  We  want 
to  show  that  the  functional 

m- 1 

p(2j-l). 


co  co  m-1 

Rf  =  f  f(x)cos  xt  dx  -  Yi  H  f(v)  -  Yj  B 

A  A  V  .  *  * 


'(0) 


v=0 


j=l 


with  the  coefficients  H  =  H^2^,  B.,,  ,  -  B[2m]  t  as  given  in  (5.42) 

v  v,  t  2j-l  2J-1,  t 


-42- 


#1183 


[or  as  given  as  the  expansion  coefficients  of  the  rational  function 
(5g  1)  where  the  are  defined  by  (5.  ?)}  has  the  property 

Rf  =  °  If  f<  S2mn  Lj(R+) 

Suppose  ft  S  n  L,(R+)  and  let 
m  l 

00 

f(x)  =  2  crQ2m(X"r) 

r=-  oo 

be  the  expansion  in  terms  of  the  forward  B- splines  of  degree  2m  -  l. 
By  Lemma  2,  we  know  that  f(x)  t  Lj(R+)  implies  that 


(5.55) 


The  partial  sums 


Y  IcJ  <  oo. 
r-0 


yx)  i  crQ2m<*-i  (k  = 


r=-oo 


have  the  additional  property  that 


(5.  56) 


*k(x)  =  0 


if  x  £  2m  +  k. 


Moreover,  f  (x)  =  f(x)  if  x?0  so  that  since  f(x)  t  n  L  (R+) 

K  2m  1  ' 

and  (5.  56)  holds,  we  conclude  that  f^(x)  t  S2m  n  Lj(R+)  for  integer 
k  =0.  Using  the  properties  of  the  functional  (5.  2)  we  obtain 


(5.  57) 


oo  oo  Hi-*  j 

/  yx)  cos  xt  dx  .  £  Hy»)  +  £  B  ^>-‘>(0). 


Observing  that  each  f^(x)  is  dominated  by  the  function 
00 

r--oo 
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which  is  sumraable  on  R+  by  Lemma  2  and  (5.  55)  we  see  that  on 
letting  k  -*  oo,  the  relation  {5.  57)  goes  over  into  the  desired  relation 

00  00  Ifr- 1 

/  f(x)cos  xt  dx  =  2)  H  f(v)+  Yi  B?i  t  f{2hl)(0). 

0  v=0  v  j=l 

This  completes  a  proof  of  Theorem  7. 


6.  The  sine  transform  (21)  for  n  =  2m.  In  §3-5,  we've  considered 
w(x)  =  cos  xt  and  n  =  2m.  We  now  want  to  consider  the  weight  function 
w(x)  =  sin  xt  and  indicate  the  modifications  in  these  previous  sections 
that  allow  us  to  prove  the  following 


Theorem  8.  Among  all  q.  f.  of  the  form 


(6.  1) 


/  f(x)sinxtdx=  J  H<2"°f(v)  +  £  B<2m)f(2i)(0)  +  Rf 
0  v=0  V»l  J=1  C)>X 


where  the  coefficients  satisfy 

I  h(2^)  |  <  K  for  fixed  t.  for  all  v  i  0  and  some  K, 
v>  c 

there  is  a  unique  q.  f.  with  the  property  of  being  exact.  Rf  =  0,  when- 

+ 

ever  f(x)  €  S2m  H  Lj(R  ).  This  unique  q.  f.  Is  given  bv 


00 


(5.  2)  /  f(x)sin  xt  dx  = 


+2m(t)  r  » 


0  W‘> 


Yj  f(v)sin  v 1 1 
v  =  l 


tyi-U£tl..<|,ZlftH2m-2i-l(1)  °°*  1 ,  {ttt).  .  Rf 

L  ,21+1  11  )f  <0,  +  Rf> 


j=0  t 


for  m  =  1,  2. 


-44- 


#1183 


For  simplicity,  we  write  Hy  =  H^2^,  =  B^2"^.  We  again 

shall  attempt  to  show  that 

I 

(6.3)  Rf  =  0  if  f  c  S.  fl  L,(R+) 

2m  l 

by  enforcing  this  requirement  for  the  sequence  of  forward  B- splines 
of  degree  2m  -  1  given  by  (3.  4).  Upon  substituting  f(x)  =  Q(x-r)  in 

I 

(6.  1)  we  have  the  sequence  of  relations 
r+2m 

(6. 4)  f  Q(x-r)sin  xt  dx  =  H  Q(-r)  +  H  Q(l-r)  +  . . , 

0  u  i 

m-1 

+  Hr«.2^10(2m-1>  +  2iB2jC3<  ”(0), 

(r  =  -2m+l,  -2m+2,:. . . ,  -2,  - 1) 

and 

r+2m 

(6.  5)  /  Q(x-r)sin  xt  dx  =  Hr+1Q(1)  +  Hr+2Q(2)  +  . . .  + 

(r=  0,  1,  2, ... ) 

which  are  the  analogues  of  (3.  8)  and  (3.  9),  respectively.  Again, 
we  use  a  generating  function  approach  and  observe  that  the  righthand 

i 

side  of  (6.  4)  and  (6.  5)  is  equal  to  the  coefficient  of  xr+2m  1  in 


oo  2m- 2  i 

(6.6)  (  X  H /)(  Yj  Q(2m-l-v)xV) 

v=0  v  =0 

2m~2  m-1 

+  E  t  Y  (2m-  l-v)Jxv. 

v=0  j=l  3 

Similar  to  our  approach  in  §4  then,  we  use  (4.  10)  and  (4.  12)  of  Theo¬ 
rem  5  to  obtain  (6,  6)  in  the  form 
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(6*7)  iEV  }  (2m~  1)!  +  ^  (2m-2j-l)I 

This  is  the  analog  of  (4.  23). 

We  consider  the  left  side  of  relations  (6.  4)  and  (6.  5),  and  define 
r+2m 

F  -  =  /  Q(x-r)sin  xt  dx  (r=-2m+l,-2m+2, . . . ,  -1,  0) 

i  0 


(6.8) 


r+2m 


r+2m-l 


=  /  Q(x-r)sin  xt  dx  (r  =  1, 2, ...). 


If  we  integrate  the  right  side  of  (6. 8)  for  (r  =  -2m+l, . . . ,  -i,  0)  by 

{ 

parts  2m  -  1  times  and  follow  the  same  procedure  we  used  in  §4,  we 

i 

get  the  following  analog  of  (4.  31) 

(6‘  «  W  i  *  ^«-r»  +^"(-r)  +-^Q,2m-2,(-r)] 

.  .m  r+2m-l  .  , 

— E  C-l)  T Y)sin(r+2m-i)t  (r=-2m+l, . . . , -1,  0) 

t  ‘  i=0 

or  by  letting  j  =  r  +  2m  -  1  and  using  (4. 10),  we  obtain 

i 

(6.  10)  (-l)mt2mF.  =  -  £  (-l)i(2|n)sin(j+l-i)t  -  tQ(2m~2)(J+l) 

3  i=0 

+  t3Q(2m_4)(j+l)  -  ...  -  (-l)mt2m"3Q!,(j+l)+(-l)mt2m_1Q(j+l) 


(j  =  0, 1, . . . ,  2m-l) . 

If  we  consider  (6.  8)  for  r  =  1,2,...  and  again  integrate  by  parts 
2m  -  1  times,  we  get  analogous  to  (4.  36)  the  relation 
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(6.11)  <-l)mt2mF  =  -  E  (-nV.Vin  (J+l-i)t 
3  1=0 


(J  =  2m,  2m+l,. . . ). 


From  (6. 10)  and  (6. 11)  we  find  that 


oo  .  f2m-l  j  .  _ 

(6.12)  -(-1)  tmX  V  =<  E  [  t  (“D  (2r)sin0+l-i)t]xJ 

j=0  3  \  j=0  1=0 

+  E  l  E  (-dV 2“)8in(J+l-i)t]A 

J=2m  1=0  j 

2ni“2  /a  .  *  2m-2  . _  . .  . 

+  t  E  Q  2m~2  0+l)x  -  t3  E  O  2m"4  0+l)xj  + ... 
j=0  i=0 

-  (-DmtZm"1  E  Qo+i)xJ. 
j=0 

We  define 

oo 

(6.13)  t(x)  =  E  [sin(v+l)t]xv 

v=0 

2tn 

and  note  that  (1-x)  t(x)  Is  precisely  the  term  In  curly  brackets  on 
the  right  side  of  (6. 12),  so  that  by  using  (4. 12)  of  Theorem  5,  we  can 
write  (6. 12)  as 


(6. 14)  -(-l)mt2m  E  VJ  =  U-*)2"1^)  +  t(l-x?m''2P  (x) 

j=0  3 

,3  ,  .  .2m- 1 

i  , c m- 4 _  #  . .  ,  .  .in  t  ^ 

-jjd-x)  P2(x)+...  -(-1)  ii^7-)|P2m_2M. 

°°  t 

Equating  (6.  7)  and  E  F  *  as  determined  from  (6.  14),  we  see  that 

J=0  3 

we  require 
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...  t.  .  42m- ,  .  .  .  .2m—  2  t^  . .  2  in-  4  . 

(6.15)  J i(x-l)  r(x)  +  ttx-l)  -  jj- P2(x)(x-1)  +  ... 


+  (-l) 


m- 1  t 


2  m —  1 


(2  m- 1)1  2  m- 2 


(x)} 


«  w  P?m,(x)  m-l  (x-l)2JP  ..  .(x) 
-  /  yTT  -irV  y  2  m- 2  _  y  q  _ 2m- 2j- 2 

’Cy""v  1  (2m- 1)!  +  ^  2j  (2m-2J-l)! 

00 

Solving  (6. 15)  for  (  Y.  H  xv)  gives  the  final  relation 

v 

v=0 


00  ,, 

(‘•>‘1  =  pkr((x- 

v=0  r2m-2  [  t 


..2m-.  . 
1)  t(x) 


m-l 

■l 

1=0 


m-l 

-  I  B2J 

j=l 


(2J+1)!  2j' 

:X-1|2|p2m-21-2(X) 

(2m-2j-l)! 


This  is  the  analog  of  (5. 1).  Our  derivation  of  (6.  16)  evidently  estab¬ 
lishes  the  following 


Proposition  2.  The  coefficients  H  =  H^2^,  B  =  of  the 

V  V,  t  J  jj  t 

most  general  functional 

00  00  m- 1  . 

(6.17)  Rf  =  /  f(x)sinxtdx-  J  H  f(v)  -  £  B  f(  J,(0) 

0  v=0  v  j=l  2) 


vanishing  for  the  function  (5.  5)  are  the  expansion  coefficients  of  the 
rational  function  (6. 1 6)  where  the  B  (1=1, . . . ,  m-l)  are  chosen 
arbitrarily. 
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1C5f'Y‘>^.<f 


We  again  wish  to  determine  the  unique  functional  (6. 17)  having 

bounded  coefficients  .  From  the  form  of  t(x)  given  by  (6. 13), 

we  note  that  t{x)  converges  for  |x|  <  1.  Then  in  a  similar  manner 

to  the  discussion  in  §5  we  observe  that  the  coefficients  H  will  be 

v 

bounded,  if  and  only  if  the  coefficients  B  (j  =  1, . . . ,  m~  1)  can 
be  chosen  to  satisfy  the  equations 


m-1  (X  -l)1 2JP  (x  )  m-1 

<6- 181  *  ■  -t^r«vi>2m;w 

m-l  (-i)j t2^1 

♦J,  (2j+l)!  VVV11  )  <v  = 

where  the  X  (v  =  1, . . . ,  m-l)  are  the  zeros  of  P,  lx)  less  than 
v  2m- 2 

one  in  absolute  value.  So  we  need  only  show  that  the  determinant 

,  ,  >-i)21p,  ,,  ,<\) 

(6-i9>  |av,'  - 

(v=l, ...,  m-lj  j  =  1 , . . . ,  m- 1 ) . 

That  this  is  the  case  is  evident  from  the  expression  (5.  40)  if  we  choose 
the  special  case  I  U  I'  =  {2,  4,  6, ... ,  2m- 2  }.  This  establishes  the 
existence  of  a  unique  functional  r-  of  the  form  (6.  17)  of  Proposition  2 
such  that  the  sequence  {H^  }  is  bounded.  The  remainder  of  a  proof 
of  the  first  part  of  Theorem  8  is  essentially  the  same  as  the  proof  of 
Theorem  7  so  we  may  omit  it. 

We  observe  that  the  functional  Rf  determined  from  (6.  2)  is  of  the 
proper  form  (6. 17)  where 
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(6-20>  no-T'1-  ^,  m 


COS  ~  4;  (t) 

Hv— ^ 


(v  =  ly  2j  !  «  .  ), 


By  (1.  3)  and  (1.8)  we  find 


(6.21) 


sin  t/2 1  2m _ 1 


<  OO  (v=l,  2,  ...  ) 


so,  for  fixed  t,  the  sequence  {H^}  as  given  by  (6.  20)  is  bounded. 
Once  we  show  that  Rf  as  given  by  (6.  2)  vanishes  for  the  sequence 
of  splines  (5.  3),  the  unlcity  established  in  the  first  part  of  the  theorem 
will  complete  the  proof  of  Theorem  8.  This  latter  task  is  accomplished 
similar  to  the  cos  xt  case  of  §5,  showing  that  Rf  =  0  for  any 
m  =  1,  2, . . .  for  (r  =  0,  1,  2, . . . )  and  looking  at  the  particular  cases 
of  m  =  1,  2  to  show  that  Rf  =  0  for  (r  =  ~2m+l,  -2m+2,  . . . ,  -1). 

We  omit  the  details. 


7.  The  even  degree  spline  case,  n  =  2m- 1.  In  §3-5  and  §6  we 
considered  w(x)  =  cos  xt  and  w(x)  =  sin  xt,  respectively,  for  the 
odd  degree  spline  case,  n  -■  2m.  Here  we  shall  consider  the  same 
weight  functions  but  take  n  =  2m-l  and  prove  the  following 


Theorem  9.  Let  m  =  2  or  3.  Among  all  q.  f.  of  the  form 


(7.  1) 


f  f(x)cos  xt  dx  =  £  H(2™-1)f(v)  +  £  B<2m“ 
0  v=0  V>  1=1 


(2m-  1)  ,(2j-  1) 


(0)  +  Rf 


(7.2)  /  f(x)sin  xt  dx  =  J  H(2™~  l)i(v)  +  £  B(2m"  1}  f(2J)(0)  +  Rf 

z  *  v  .  t  t  ,  L 1 ,  t 
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_ 


where  the  coefficients  satisfy 


|H(2m-l) 


.  I  <  K  for  fixed  t,  for  all  v  10  and  some 

t 1 


there  is  a 


with  the  property  of  being  exactf  Rf  =  0,  when¬ 


ever  f(x)  c  j  n  L^(R+).  This  unique  q.  f.  is  given 


oo  >\f  (t)  T  co 

(7.  3)  /  f(x)cos  xt  dx  ,  -M(0)  +  Yj  f(v)cos  vt 

0  92n>-r7  v  =  l 


n£l  /  ,J  i(t)  r>\  ^\ 

fj  (0)+  Rf 

J  =  1  t  J  92m-l'  7 


00  l}j  (t)  I  00 

(7.4)  / f(x)sin  xt  dx  =  --  1 :  .  i  Y  f(v)sin  vt 

0  92m-r  7  I  v  =  l 


/  l  J  j(i)  COS 

+  E  -^tl—  1  ,  - -]f(  ”(0)  +Rf. 

J=0t2,+1  *2m-l<t) 


We  shall  indicate  the  modifications  in  §3-5  necessary  to  prove 
this  theorem  for  the  weight  function  cos  xt.  The  case  of  sin  xt  then 
follows  as  §6  did  for  the  case  n  =  2m.  For  simplicity,  we  write 

H  i  =  H^2™  B  =  b!2™  ^ .  We  again  attempt  to  show  for 

V  l  uj”  I  tj"  lj  v 


(7.  l)  that 
(7.  5) 


Rf  =  0  if  f  €  S_  n  L,(F  ) 
2m- 1  1  ' 


by  enforcing  this  requirement  for  a  sequence  of  B- splines.  This  time 


we  choose  to  use  the  sequence  of  central  B- splines  of  degree  2m  -  2 
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(7.6)  (M2m_1(x-m-r)}  =  (M(x-m-r)}  (r=-2m+l,~2m+2, . . . ) 

whose  knots  are  at  the  points  (v  +*~).  v  an  integer.  We  write  these 

& 

B- splines  in  the  form  (7.  6)  to  make  the  analogy  with  §3-5  clearer. 

By  (1.  5)  we  have  the  explicit  expression 


.  2m- 1 

*7,7)  M2m-l(x)  =  (2m-2)l 


v  ,  ..i  . 2m-l . ,  .  2m-l 
L  (-D  (  <  )(x  +  — r~  -i) 


2m-2 

+ 


tjpon  substituting  f{x)  =  M(x-m-r)  in  (7. 1)  and  noting  the  requirement 


(7.  5),  we  obtain  the  sequence  of  relations 


r+2m-~ 

(7.8)  /  M{x-m-r)cos  xt  dx  =  H0M(-m-r)  + Hj  M(l-m-r)  +  . . . 


*  Hr+2m-l^m-^  1  m— r) 

(r  =  -2m+l, . . .  ,-l) 

and  j 

r+2m--2 

(7.9)  f  l  M(x-m-r)cos  xt  dx  =  H +1M(l-m)+H +2M(2-m)  +  . . . 

r+i 

+  ^r^2m- ^  =  0»  1>  2, . . . ) 

which  are  the  analogs  of  (3.  8)  and  (3.  9),  respectively.  We  again  em¬ 
ploy  a  generating  function  approach  and  note  that  the  right  side  of 
(7.  8)  and  (7.  9)  is  equal  to  the  coefficient  of  x  in 

oo  2m-2  2m-2m-l 

(7.  10)  (  Ys  H  xV){  Z  M(m-l-v)xV)  +  £  [  J  B 

v=0  V  v=0  v=0  J=0  J 

In  order  to  simplify  the  two  polynomials  in  (7.  10)  we  note  that 
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(7,11)  M(k)(x)  =  (-l)k  M(k)(-x) 


as  can  be  verified  from  (7.  7).  With  this  substitution  and  (4.  14)  of 
Theorem  5,  we  obtain  (7,  10)  in  the  form 

oo  T  (x)  m—  1 

n  I  V  TT  /l4l2ia“2  2m~2  V  n  ,i.2m-2J-l 

(7.12)  (  IjHx)(2)  i7fn_7U  “  L*7i_t(7) 


(2m-2)! 


y  _  \  2m-2j-l 

■  AViV 


T2m-2,-i<*> 


(2m-2j- 1)! 


This  is  the  analog  of  (4.  23). 

We  consider  the  left  side  of  relations  (7.  8)  and  (7.  9),  and  define 


(7.13) 


r+2m- 


r+2m- 


r+2m~— 

=  /  M(x-m-r)  cos  xt  dx 

0  (r  =  -2m+l, -2m+2, . . . , -2, -1) 

r+2m-“ 

=  f  L 

1  r.l  M(x-m-r)cos  xt  dx  (rs  0.  1.2....). 


(r  :  0j  lj  2^ , , , ) , 


If  we  integrate  the  right  side  of  (7.  13)  for  (r  =  ~2nvflt  -2m+2,  ,  -1) 

by  parts  2m- 2  times  and  follow  the  same  procedure  we  used  in  §4,  we 
get  the  following  analog  of  (4,  31) 


(7.  14)  Fr+,m_  j  =  -  [“  M* (-m-r)  -  —  M"  (-m-r)  +  . , . 


,2m- 2 


..(2m-3).  .. 

M  (-m-r)J 


.m-2  r+2m-l 


“  Yj  (-l)1(2l^i)sin(2m  +  r-  1  +~-  i)t 


(r  =  -2m+l, . . . 1) 
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or  by  letting  j  =  r  +  2m  -  1  and  using  (7, 1 1)  we  obtain 


(7.15)  (-i)n"1t2m~1F  =  S  (-l)1!2"”1) sin(j+^-i)t- 

3  i=0  1  c 


tM(2m-3)(J+i_m) 


+  t3M(2m“5)(j+l-m)-  ...  +  (-l)m“1t2m”3M,{j+l-m) 


()  =  0, 1, . . . ,  2m-2). 

If  we  consider  (7.  13)  for  (r  =  0,  l,  2, . . . )  and  again  integrate  by 
parts  2m- 2  times,  we  get  analogous  to  (4.  36)  the  relation 

(7.16)  (-l)lT>"1t2m“1F  =  Z  (-l)i(2?”1)sin  0+4-  i)t 

3  i=0 

(J  =  2m- 1,  2m, . . . ) . 

From  (7.  15)  and  (7.  l6),  we  find  that 


(7.  17)  (- 1)*”“ 1  t2m~ 1  £  = 

J=0  3 


=  < 


( 2m-2  J  ,  ,  ,  ,  . 

Z  [  t  (-D^2?  )sin(j+j-i)t]xJ 

j=0  i=0  c 


+  Z  ['ZM^^^sinO^-iM 

j=2m- 1  i=0  J 

2m-2  .  .  ,  2m-3  ..  . 

t  Z  M  ‘n  (J+l-m)x^  +  t  Z  M  m  (j+l-m)x  - 
1=0  1=0 


2  m- 2 

nr1^3  z 

1= 0 


M*(J+l-m)x*. 


We  define 
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(7.18)  U(x)  =  X  (sin  (v  +  4)t]xV 

v-0 
2  1 

and  note  that  (1-x)  ^(x)  is  precisely  the  term  in  curly  brackets 

on  the  right  side  of  (7.  17),  so  that  by  using  (4.  14)  of  Theorem  5,  we 
can  write  (7.  17)  as 

oo  fl-x)2m“3T  fx) 

.  .,m-K2m-l  v  „  1  .2m-l . t.'1  V  ' 

(7.19)  (-1)  t  A  F  x  =  (1-x)  U(x)-(t) - — - 

j=0  J  2  1! 

(l-x)2m"5T  (x) 

+  <I>  - 3!  ’  "• 


,  .m-1  t  2m-3  (I_x)T2m-3(x) 

+  <_1)  (2!  (2m-3)!  ' 

oo 

Equating  (7.  10)  and  F  x^  as  determined  from  (7. 19)  and  then  solv- 
oo  J=0  J 

ing  for  ^  Hxl  gives  the  final  relation 


00  ,...,2m-2  f,  ,  4m-l  ■»  i 

<7-*» 

v=0  2m-2  t 


m-1  ...  .  (x-1)2J-1T  fx)  'I 

V  P  /i/m-2j-l  1  2m- 21-1  | 

“j"j  2J-1' 2  *  (2  m— 2  j—  1 ) !  j 


This  is  the  analog  of  (5.  1),  Our  derivation  of  (7.  20)  evidently 
establishes 


Proposition  3.  The  coefficients  =  H 


(2m- 1)  _  (2m- 1) 

v,  t  >  2j-l  "  2j-l,t 
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of  the  most  general  functional 

oo  oo  m-1  . 

(7.21)  Rf  =  /  f(x)cos  xt  dx  -  SHf(v)-S  B  f*  *  '(0) 

0  v=0  V  J=1  2J"1 

vanishing  for  the  functlons(7.  6)  are  the  expansion  coefficients  of  the 
rational  function  (7.  20)_wherejhe  B  (J  =  1, . . . ,  m- 1)  are  chosen 
arbitrarily. 

We  again  wish  to  determine  the  unique  functional  (7.21)  having 

bounded  coefficients  H^.  From  the  form  of  U(x)  given  by  (7. 18),  we 

see  that  U(x)  converges  for  |x|  <  1.  So  that  Just  as  in  §5  we  note 

that  the  coefficients  H  will  be  bounded  if  and  only  if  the  coefficients 

v 

®2J  1  M  ♦  •  •  >  m_1)  can  be  chosen  to  satisfy 


(7.  22) 


m~  1 


.J  2m-2J-  1  ^  T2m-21-l(V 

'  2  ’  (2m-2j-  1)  ! 


.  ..m-1  ,  ,  m-1  .  .  .  T  _(\  )(\-l) 

”<V  + £/-»’<!>  -  (Zi-in 


2m-2j-l 


(v  =  1,  2, . . . ,  m-1) 

where  the  \  (v  =  l,...m-l)  are  the  zeros  of  T,  _(x)  less  than 

one  in  absolute  value.  Lemma  8  of  [9,  p.  182  J  had  guaranteed  that 

T_  _(x)  was  a  reciprocal  polynomial  which  had  only  simple  and 
2m- 2 

negative  zeros  X.  ,  x.^, . . . ,  \  ^  that  we  may  label  to  satisfy  the 

conditions  (5.  5)  and  (5.  6).  So  we  need  only  show  that  the  deter¬ 
minant 
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1  1  (^-D2  1  T,m  „  I<k  > 

(7  23)  Ia  I  =  |(l)2mr2J“1  — K - -ImilhJ—iL  {  J  0  . 

(7. 23)  vj  ‘  '  2  *  (2m-2j-l)l  1 

We  haven’t  proved  a  ler.ima  similar  to  Lemma  4  of  §5,  but  note  that 
the  determinant  in  (7.23)  for  the  cases  m  =  2  and  m  -  3  takes 
the  forms 


(7.  24) 


~<X-1)(X+1) 


(7.25)  ^(X1-1)(X2-1)(X1+1)(X2+1)(X1-X2)(1-X1X2) 

respectively,  and  so  by  (5,  5)  the  condition  {7.  23)  is  satisfied. 

This  establishes  the  existence  of  a  unique  functional  Rf  of  the  form 
(7.  21)  of  Proposition  3  such  that  the  sequence  (H^  }  is  bounded.  The 
remainder  of  the  proof  of  the  first  part  of  Theorem  9  is  essentially  the 
same  as  the  proof  of  Theorem  7  and  we  omit  it. 

We  observe  that  the  functional  Rf  determined  from  (7.  3)  is  of  the 


appropriate  form  (7.  21)  where 

!7-26) 


♦2m-l(t> 


cos  vt  (v  =  1,  2, . . . ) 


and  by  (1.  3)  and  (l.  8)  we  again  have  this  sequence  }  bounded. 
It  is  a  straightforward  procedure  to  show  for  m  =  2,  m  =  3  that  Rf 
as  given  by  (7,  3)  vanishes  for  the  sequence  of  splines  (7,  6),  so  that 
the  unicity  established  in  the  first  part  of  Theorem  9  completes  the 
proof  of  the  theorem. 
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8.  Explicit  forms  for  the  q.  f.  of  Theorem  2.  Theorem  2  is  es¬ 
tablished  similar  to  the  way  the  first  parts  of  Theorem  7  and  8  were 
proved.  The  only  change  needed  is  to  take  I  U  I'  =  {1,2 , , . . ,  m- 1 } 
instead  of  what  we  used  before.  For  instance,  for  the  cos  xt 
case,  we  use  the  same  left  sides  of  (3.8)  and  (3.9),  but  we  have  to 
modify  the  right  sides.  Theorem  5  readily  allows  us  to  do  this  and  we 
get  a  similar  expression  to  (4. 42)  but  now  in  the  coefficients 
B ..  B,t , . . ,  B  ,  .  Lemma  4  for  the  choice  I  U  I'  =  {l,  2. . . . ,  m- 1  } 
enables  us  to  establish  Theorem  2  just  as  the  first  part  of  Theorem  7 
was  proved. 

For  the  weight  function  cos  xt,  the  q.  f.  (5. 42)  gives  explicit 
expressions  for  our  present  cases  m  -  1  and  m  =  2.  We  want  to  get 
a  q.  f.  similar  to  (5.42)  for  m  =  3,  the  quintic  spline  case,  in  a  form 
particularly  amenable  to  computation.  We  shall  find  that  the  form  we 
do  obtain  is  precisely  (5. 42)  with  the  exception  that  f"'(0)  is  re¬ 
placed  by  S'"(0),  the  third  derivative  of  a  particular  interpolating 
spline  to  f(x),  evaluated  at  0.  The  expression  for  S"'(0)  involves 
the  values  f"(0),  f'(Q),  f(0),  f(l),  but  not  f'"(0). 

For  the  sin  xt  case,  the  q,  f.  (6,  2)  gives  our  desired  q.  f.  when 
m  =  1.  For  m  =  2,  the  only  change  we  make  in  (6.  2)  is  to  replace 
£"{G)  by  S"(0),  the  second  derivative  of  a  particular  spline  interpolant 
to  f(x),  evaluated  at  0.  Here  S"(0)  is  expressed  in  terms  of  the 
values  V  (0),  f{0),  f(l),  ...  but  not  f  "(0).  There  is  a  similar  q.  f. 
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J  wwwww-w **++—*,, 


(iv)  (iv) 

when  m  =  3  that  we  get  from  (6*  2)  by  replacing  f'  (0)  by  S  (0). 
We  summarize  the  foregoing  in  the  following 


Theorem  Z\  Suppose  f«  Lj(R+).  Let.  S(x)  c  Lj(R+)  be  the 
unique  spline  of  degree  2m  -  1  for  x  =  0  with  knots  at  x  =  1,  2,  3, . 
satisfying  the  conditions 

(8.1)  S(v)  =  f(v)  (v  s  0,1,2,...) 

(8.2)  S(J,(0)  =  f0)(0)  (J  =  l,2,...,m-l). 

i 

For  m  =  1,  2,  3,  the  unique  a,  f.  of  Theorem  2  are  given  explicitly  by 


00  4l  (t)  f  00  1 

(8.3)  /  f(x)cos  xt  dx  =  — Jrf(0)  +  2j  f(v)cos  vt\ 

0  *2m(t'  \ 2  v  =  l  J 


and 


[m/2]  ,  j  ^.(M-  ,,(*)  m  n 

+  l  ^n--1 .  ;( 1  <°) 

(  =  1  t  '  +2m(t) 


m-  1 

l 


J=(m/2]+l  t 


2J 


j 

(1- 


.(t) 


(-l)^f ,  _  W'^zm-zr'  ,  c(2j~  1) 


*2m(l) 


3  s' 


(0)  +  Rf 


00 


(t)  f  00 

(8.4)  /  f(x)sin  xt  dx  = ^  £  f(v)sin  vt^. 

0 


*2m(t) 


v=l 


J 


l^r) 


,(t)COS- 


2  .  J  .(‘/wo  .... 

,  V  izJJL  f ,  -JJ _ 2m- 2]-  1 - Z_  l  f'2J).Q% 

+  b  2j  +1  1  ^  ft)  1  1  ' 


j=0  t 
m-1 


*2m(t) 


izii.fi  .  _JJ - Z-7~ ■ - - ] Sf2j>(0)  +  F 


2J« 


*2m(t) 


I 


We  also  state  here  for  reference  in  §10  the  following 

I 

Corollary  l.  Let  f(x)  and  S{x)  be  as  given  in  Theorem  Z',  Then 
the  Qi  f. 


/  f(x)sltxdx  =  ^2^  I  i((0)  +  J  f(v)ett” 

o  92m' 1  \  W  =  1 


'j«l  t2J  ^2m 


(t) 


Jf^CO) 


n»-i  ,  n J  '  ?4(t)  (?i  n 

+  ^  [1-— J— 2-~-2i-]S(2J  U(o) 


i=[m/2]+l  t2J 


+  i< 


r[(nrf)/2]  J  „  ,(t)cost 


r  tl  -  1 t<2»(0, 


1=0  t 


♦2ra(t) 


m-1 


Vi  /  (t)cos- 

X  - '•i8  <°>y+R1 


)=[(m-l)/2]+l  t2)+1 


Is  the  unique  q.  f.  exact  whenever  f(x)  c  3,  n  LfR  ). 

2  m  1 

We  discuss  in  detail  the  cubic  case,  m  =  2,  of  the  q.  f.  (8. 4), 
that  is,  the  q.  f, 

00  vji4(t)  »  4»^(t)cos  t/2 

(8.5)  f  f(x)sin  xt  dx  =  7-7^  £  *(v)sin  vt  +  ~[ l  - 
0 


*4(t)  v=l 


<M*> 


]f(0) 


.1  r.  "1 


4>,(t)  cos  t/2 


[1- 


*4(t) 


]S"(0)  Rf  . 
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t 

* 

f 

i 

<i 


i 

i 


j 

v 

i, 


How  to  extend  the  method  to  m  =  3,  and  in  fact  larger  m,  and  (8.  3)  will 
be  clear.  We  observe  that  by  the  second  part  of  Theorem  8,  the  exact¬ 
ness  of  (8.4)  for  f(x)  <  S2m  t;  L1(F+)  is  already  established  since 
for  such  f  we  have  S”(0)  =  f”(0).  What  remains  to  show  is  that 

i  O  __  _  \ 

the  coefficients  H'  ,  of  the  q.  f.  (17)  of  the  Introduction  are  bounded. 

v,  t 

An  important  point  is  that  we  do  not  want  to  use  the  ’'natural1' 
semi-cardinal  cubic  spline  interpolant.  This  is  of  course  best  in  the 
sense  that  it  minimizes 

00 

/  [F"(x)J  dx 
0 

among  all  functions  F(x)  that  interpolate  f(x)  at  x  =  0,  1,  2, . . .  but 
it  is  not  a  good  approximation.  Rather  we  use  the  "complete"  semi¬ 
cardinal  spline  approximation  where  also  f ' (0 )  is  assumed  known 
and  is  matched  by  the  cubic  spline,  that  is,  (8.  2)  holds. 

We  note  that  Lemma  4  guarantees  the  unicity  of  the  interpolating 
spline.  The  interpolating  spline  S(x)  is  given  by  the  spline  interpo¬ 
lation  formula 

00 

(8.  6)  S(x)  =  Yj  f(v>L  <x>  +  f'  (0)  A(x) 

-  V 

v  -0 

where  the  fundamental  functions  L  (x)  and  A(x)  satisfy 

v 

L  (v)  =  1 
v 

(8.7)  Lv(h-)  =0  if  v/u  (v  =  0,  1,  2,  .  .  .  ) 

L'(0)  =  0 
v 

and 
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(8.  8)  <v  =  0,  1,  2, . . . ) 

A'(0)  =  1 

respectively.  In  order  to  construct  Ux)  and  A(x)  we  shall  use 
two  other  important  cubic  splines. 

One  is  the  fundamental  function  L(x)  of  cardinal  cubic  spline 
interpolation,  i.  e. ,  L(x)  is  a  cardinal  cubic  spline  satisfying 

{1  if  v  =  0 

0  if  v/0. 

In  terms  of  the  cubic  B-  spline  M(x)  =  M^fx),  we  have  explicitly 

0°  I  .  I 

(8.  10)  L(x)  =  n/3  £  X  1,1  M(x~j) 

J=-oo 

where 

(8.  11)  X  =  -  2  +  =  **  .  268  . . . 

is  the  root  of  least  absolute  value  of 

P2(\)  =  X2  +  4X+  1  =  0  . 

To  prove  (8.  10)  we  note  that,  for  v  =  1 

L(v)  =  n/3  Yj  X^M(v-J)  =  ^3"(XV+1  +  4X%  XV_1) 
i-v- 1 

=  "/I  •  “  ’  XV-1(X2  +  4X  +  1)  a  ) 

and 

L(0)  =  n/T  .  •  (X+4+X)  =  “  (2+X)  =  1 

so  that  the  unicity  of  a  bounded  L(x)  satisfying  *8.  9}  implies  that 


(8.  10)  is  correct. 
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The  other  cubic  spline  we  use  is  the  decreasing  cubic  eigenspline 


(8*  12) 


. 

S  (x)  =  Yi  *  M(x-j) 
j=-00 


which,  from  §5,  is  a  cardinal  spline  satisiying 


(8.  13) 


(8.  14) 


Now  if  we  write 


(8.  15) 


Sj(x)  =  o(Mx) 


Sj(v)  =  0 


A(x)  =  Sj*  (01  S1(X) 


as  x  -*■  oo 


for  all  v . 


and  use  (8.  14),  we  see  that  (8.8)  is  satisfied.  From  (8.  12), 

(8.16)  Sj'fO)  =  s  \jM'H)  rjfX-X-1)  =  s fl 

because  M'(-l)=-~,  M*(0)  =  0,  M'(l)=-4and  M’  (j)  =  0  for  all 

u  z 

other  j.  From  (8.  15)  and  (8.  12) 


(8.  17) 
Writing 
(8.  18) 
we  see  that 


1  00  1 

AW  =  j?  E  XJM(x-J)  for  X  £  0. 
J=-l 

00 

A(x)  =  Yj  Y,M(x-j) 

J=-l  J 


(8.  19) 


V  =  T=XJ 
Vj  >f3X 


We  shall  use  the  notation 


(8.  20) 


We  ha'fe 


(8.  21) 


Lv(x)  = 


00 

l  «, 

j=-l  U 


M(x-j) 


L0(x)  =  L(x) 


(j  =  -1,  0,  1, ...  ). 


for  x  =  0 . 


for  x  =  0 


since  by  (8.  8)  and  1/(0)  =  0,  (8.  7)  is  satisfied.  Using  (8.  10)  and 
(8.  21),  we  find 


(8.  22) 

It  is  easily  verified  that 


c  =  v 3 

Cj,0 


(j  =  -1,  0,1,...). 


(8 


.  23)  Ux)  =  L(x-v)  -  Sj(x)  (v=l,2,...) 


satisfies  the  conditions  (8.  7).  From  (8.  10)  and  (8.  16) 

/mv‘fl.  'v“1‘ 

1'"'  X-X* 

and  therefore  by  (8.  23),  (8.  10)  and  (8.  12) 


.  KIsL  ..  d x*"1)  .  j,  (v  .  ,  - 

S«(0)  -  .  -1  -  ^3  \  (v  _  1, 2, ... ) 


L  (x)  =  ^3  E  X  ^  M(x-v-j)  +  n/3  XV  ^XJM(x-j) 

V  J  J 


—  hi  —  00 

=  *^3  E  X,J“V,M(x-j)  +  V 3  XV  X^M(x-J)  =  E  M(x-j) 


i,v 


where 
(8.  24) 


J=-l 

for  x  £  0 


CJ,v  =  ^3(X*j"V^  +XJ+V)  (v  £  1,  J  =  -1,0,  1,...). 


We  can  now  use  the  interpolating  spline  S(x)  given  in  (8.  6)  to 
determine  from  the  q.  f.  (8.  5)  what  the  form  of  the  coefficients 
(4) 

IT  t  of  q.  f.  (17)  is.  Differentiating  each  side  of  (8.  6)  twice 

and  substituting  0  for  x  gives 

00 

(8.25)  S"(0)  »  J  f(v)L"(0)  +  f’(0)A"(0) 

v-0  V 


so  that 
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.  «|»  (t)cos-j 

(8- 26)  Hv  =  ITrn sln  vt  jt1  3  Lv (0)  (v=  *>  2>  •  ♦  •  > 

4  t  m  ’ 

From  (8.  20)  we  have 

00 

(8.27)  L  "(0)  =  Yj  c,  M”H)  =  c  ,  -2c  +  c, 

v  jifj  J,v  -],v  0,v  l,v 

since  M"(- 1)  =  M"(l)  =  +1,  Mf,(0)  = -2  and  M"(j)  =  0  for  all  other 

J.  From  (8.24)  and  (8.  27),  then,  after  simplification  we  get 

(8.28)  L”(0)  =  -12V3  \  (v  =  1,2,...). 

Therefore,  for  fixed  t,  by  (6.21),  (8.  11),  (8.26)  and  (8.28)  there 
exists  a  constant  K  such  that 

|H  |  <K  for  all  v  5  0. 

v 

By  the  unicity  established  in  Theorem  2,  Theorem  2'  is  established 
for  m  =  2. 

We  also  note  that  in  terms  of  the  functions  L  (x)  and  A(x) 

v 

Just  defined,  we  have  the  following 
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f°°  .  »}/  Wcos^ 

/  A(x)  sin  xt  dx  =  -  -*-[l  -  — 1  Z  ] 

0  t  ♦4*t) 

,w  i+iW  1  «Mt)cos- 

/  L0(x)cosxt  dx=i— ,  /  L0(x)slnxtdx=jl-— »)-2]. 


-66- 


#1183 


III.  THE  LAPLACE  TRANSFORM 


We  follow  the  same  kind  of  generating  function  approach  we  used 
earlier,  the  only  modifications  coming  from  the  enlarged  class  of  func¬ 
tions  for  which  we  can  find  transforms  now  that  the  weight  function  is 
e“xp,  p  >  0. 


9.  Proof  of  Theorem  3  of  the  Introduction.  We  first  define,  for 

each  y  =  0,  the  class  of  functions 

F  =  {F(x)  :  F{x)  «  C(R+)  and  F(x)  =  0(xY)  as  x  -  +  oo  } 
Y 

and  the  class  of  sequences 


Y  =  (v  =  {yvCr  yv  =  0(vY)  as  v  "  00 ) 


We  note  that 


(9.1)  S(x)  s  S2m  H  F^  for  some  y  =  °  implies  {S(v) }  c 

so  that  by  (24)  the  functional  Rf  given  by  (23)  is  well-defined  if 

f(x)  c  S_  OF.  We  also  need  to  know  just  when  a  cardinal  spline 
2m  y 

function  of  degree  2m  -  1  is  in  F  ,  The  answer  is  given  by 

Y 


Lernma  5.  If  S(x)  «  S„  and 


(9.2) 


s(x)  =  Ecv  M2m(x-v) 


(9.  3) 


S(x)  c  F 


if  and  only  if 


(9.4) 
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{c  }  €  Y  . 

V  Y 


Proof.  First,  we  assume  (9. 4),  so  there  exist  constants  A,  N 


such  that 

(9.5)  jc  I  <AvY  for  v>N. 

v 

Let  n  >  m  +  N.  From 


*  n+m 

|S(x)|  £  £1°  |m_  (x-v)  (0)  V  jc  |  if  n  =  x  =  n+1 

v  2m  2m  v 

v  v=n-m 

and  (9.4),  we  find 

n+m 

(9.6)  |S(x)|  S  M2m«»"  A  Tj  vy£[M2  (0)*  A*  (2m+l)K]nY  ^  k|x|y 

!  v=n-m  m 

if  n  £  x  3?  n+ 1 

where  K  represents  the  quantity  in  square  brackets  in  (9.  6)  and  does 
not  depend  on  n.  So  for  any  x  >  n  >  tn+N 


|S(x)j  iK|xJY 

and  (9.  3)  holds. 

Now,  we  assume  (9.3)  and  are  to  derive  (9.4),  We  adapt  the  proof 
of  Theorem  4  of  [  1 1,  p.  18- 19]  to  our  particular  situation,,  We  observe 
first  that  in  Theorem  5  of  [  11,  p.  7]  Schoenberg  explicitly  expresses 
the  c  of  (9.  2)  in  the  form 


(9.7) 


"v1  /  J  (m)  „(2r) 

C  =  1  (-D  Yor  S  V) 


r=0 


(m) 


where  the  }  is  a  sequence  of  rational  numbers  generated  by  the 


expression 
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i' 


r 


r _  u  ,2  m  ^  (m)  2r 

‘  2  sin  u/2  '  =  v2r  U 

i"=0 


By  (9.  7)  then  in  order  to  prove  (9.  4)  it  is  sufficient  to  show  that 
«•«>  {S(2r)(v)}<  Y  (r  =  0, 1, . . . ,  m- 1) . 


By  (9.3 )  there  exist  constants  C,  D  such  that 

(9.9)  |S(x)|  =  CxY  for  x  >  D  . 

Let  R(x)  be  a  polynomial  of  degree  k  in  the  interval  [0,  l].  By 
Markov's  theorem  we  obtain  the  string  of  inequalities 

max  | R*  |  ^  2k  max  |r| 
max  |  R"|  s  2(k-l)2  max  |  R*  | 


»  ♦ 

t 


max|R(t)|  5  2{k-t+l)2  max|R(t"  1J|  (t  ^  k, 

and  putting  them  together  we  obtain 

(9.10)  max|R(t)|  5  A(k,  t)  max  |r|  (t  =  1,  ...,k) 

where  A(k,  t)  =  2t[k(k-l)  ...  (k-t+1)]2  .  Applying  (9.9)  with 
k  =  2m-  1,  t  =  2r,  to  each  of  the  polynomial  components  of  S(x)  in 
each  of  the  successive  intervals  [v,  v+l]  for  v  >  B  +  1,  we  con¬ 
clude  that 


|  S{2r)(v)  |  =  max  !s{2r)(x)|  ~  A(2m-l,2r)  max  |S(x)i 
[v-1,  v]  [v-l,v] 


=  A(2m-1,  2r).  C  max  xY  =  A(2m-  1,  2r)C  •  vV 
[v-  1,  v] 

so  that  (9.  8)  and  the  lemma  are  proved. 

#1183 


-69- 


We  now  turn  to  the  proof  of  Theorem  3  and  again  shall  indicate 


the  modifications  in  a  previous  proof,  this  time  Theorem  2,  that  leads 
us  to  our  desired  conclusion.  We  again  require  exactness  for  the 
B- splines  (3.4)  and  are  led  to  the  sequence  of  relations 
r+2m 

(9.  11)  /  Q(x-r)e~xpdx  =  HQQ(-r)  +  HjQd-rH. . .  +Hr+2m_1Q(2m-l) 

1 

+  Z  BO  {-r)  (r  =  -2m+l,-2m+2,...,  l) 

J=1  J 

and 

r+2m 

(9.  12)  /  Q(x-r)e~Xpdx  =  H^QU)  +  Hr+2Q(2)  +  . . . 

+  Hr+2m_iQ(2m-l)  (r  =  0,  1,2,...). 

The  righthand  sides  of  (9.  11)  and  (9.  12)  are  the  same  as  in  the  proof 
of  Theorem  2,  but  the  left  sides  now  have  the  weight  function  e  Xp . 
Following  the  same  procedure  as  before,  we  arrive  at  the  following 
analog  of  (5.  l) 


(9.  13)  X  H 


l  H  x%Jl!!tliL  IO_ 

v=0  v  P2m-2W  p2m 


_  -p  2m- 1  J 

[<x-U2m-^r+  l  jj 

1-e  px  J=0  J 


•  Pj^jtxJfx-l) 


J  '  j— j  J  (2m-l-J)( 


Let  R  (x)  denote  the  right  side  of  (9.  13),  where  the  B. 

2m  j 

(j  -  1,2,...,  m- 1 )  are  as  yet  undetermined.  We  recall  that  P  ^(x) 
has  the  simple  zeros  \  ,  X.  , . . . ,  ,  satisfying  (5.  5)  and  (5.  6). 

x.  Ct  01”*'  c» 
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are  inside  the 


Observing  that  for  R,  (x)  the  poles  x,,...,X 

l  m-l 

unit  circle  while  ep  and  X  .....  X^  „  are  outside,  in  view  of 

nr  *  2m- 2 

(5.  5)  and  (9. 13),  we  note  that  the  coefficients  H  will  satisfy  the 

v 

condition  (24)  if  and  only  if  the  coefficients  (j  =  1,  2, . . . ,  m-1) 

can  be  chosen  so  that  the  m-1  poles  X., ....  X  .  have  vanishing 

r  *  m-1 

residues.  By  (9. 13)  this  will  occur  if  and  only  if  the  satisfy  the 
equations 

m-1  (X-l)jP.  ,  .(X  )  .  ,  -p 

.....  V'  «  v  2m-2-J  v  1  ,.2m  e 

(9.14)  2j  B  — ;-;-~Tn7 —  =  - ~ 

1=1  1  (2m-  1-J>!  ”  i-e-pK 

V 

2m-l  j  2ii 

+  L  *fr  P2J-l(Xv,<Xv_1)  ')  (v=l,2,...,m-l). 

The  determinant  of  the  system  (9.  14)  though  is  not  zero  as  is  evident 

from  the  expression  (5.40)  if  we  choose  the  special  case 

I  U  I'  =  {l,  2, . . . ,  m-1  }.  This  establishes  the  existence  of  a  unique 

functional  Rf  defined  by  (2  3)  whose  coefficients  satisfy  (24)  and 

which  vanishes  for  the  functions  (5.  3). 

The  remainder  of  a  proof  of  Theorem  3  follows  the  same  procedure 

as  the  proof  of  Theorem  7,  where  now  we  use  F  and  e  Xp  instead 

Y 

of  L^F*)  and  cos  xt  respectively,  and  so  we  omit  it. 


10.  An 


version  of  the  a.  f.  of  Tneorem  3.  In  terms  of  the 


central  B- spline  M  (x)  of  degree  n-1,  we  define 

n  ’ 
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(10.1) 


00 


/  Mn(x)e"xpdx  =  i  (p) 
-eo 


which  upon  evaluation  gives  the  relation 


(10.2) 

We  also  define 


7  ,  ,  r  2  sinh  p/2,n 
*n(p)  =  1 - T^J  * 


(10.  3) 


*nM  = 


l  M^vle-O  =  S  Mn(v)e"vp 

V=-00  I V  I 


and  note  that  4>n(p)  Is  positive.  From  (10.  3)  we  find 


i2(p) 

<fr3(p) 

♦4<P) 

i5(p) 

^6(p) 


1 

3  4-  cosh  p 
4 

2  +  cosh  p 
3 

115  76  cosh  p  +  cosh  2p 

192 

33  +  26  cosh  p  +■  cosh  2p 
60 


and  observe  that  4>n(p)  has  the  form  of  4>n(p)  given  by  (1.7)  if  we 
replace  the  cosine  function  by  the  cosh  function.  Similarly,  if  we 
replace  the  sine  function  in  the  expression  of  4>n(p)  in  (1.  3)  by  the 
sinh  function,  we  get  precisely  4*n(p)  as  given  by  (10.  2). 

Now  we  can  state  the  following 


Theorem  3*.  Suppose  f(x)  <  for  some  y  ~0.  Let  S(x)  <  F 
be  the  unique  spline  of  degree  2m- 1  for  x  ~0  with  knots  at 
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'QStWF?  <w<5r~*w.a „ 


x  =  1,2,3,...  satisfying  the  conditions 


(10.  4) 
(10.  5) 


S(v)  =  f(v) 


(v  =0,1,2,...) 


S(J)(0)  =  f(J)(0)  0  =  1,  2, . . . ,  m-1). 


For  m  =  1,  2,  3  the  unique  q.  f.  of  Theorem  3  Is  given  explicitly  by 


°°  vn  4»?m(p)  , 
(10.6)  /  f(x)e  pdx  =  — — rl 

0 


00 


(~f(0)  +  y,  e"vpf(v)] 

V  =  1 


+  1  {1  »2](p)*2m-2l(v)]1(2)-l) 

1=1  P2' 


(0) 


«,2m(p! 


Ml 

♦  l 


m"1  in-tfiVigifri),,, 


M7M p 

m 


tli  , 

+  S  2J+1  ^  ’ 

j=0  p 


»21(p)^2m-2]-l(p)  COShl ,  ,(2J) 


_ -2J- 

*2m(p) 


K  ,#(0> 


m-1 

+  £  2J+1 

p 


-  P 

i  «l>?i(p)^?m  ?1  . (p)cosh ~  .... 

~  fl  -  - ZULU  1 - — ]  S  2  (0)  +  Rf 


*2m<p> 


where  we’ve  written  4>q(p)  =  1  for  notational  convenience. 

We  first  remark  that  as  a  result  of  Lemma  2  of  [  11,  p.  12]  we 
have  the  following: 


Every  S(x)  t  D  F  may  be  uniquely  represented  in  the  form 


(10.  ?) 


m-1 

S{x)  =  Yj 


k  =  l 
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for  appropriate  values  of  the  coefficients  a^  . 

Then,  by  (5.  22)  each  such  S(x)  also  satisfies 

S(x) «  S®m  D  Lj(R+)  for  some  r  =  0,  1, , . . ,  2m~l 
and  Lemma  4  applies.  This  guarantees  the  unlcityof  the  interpolat- 

I 

ing  spline  of  Theorem  3‘. 

We  could  follow  the  same  type  of  procedure  we  used  for  Theorems 
8  and  10  to  show:that  (10.  6)  actually  is  the  unique  q.  f.  of  Theorem  3, 
but  we  shall  not.  Instead,  we  note  that  (10.  6)  follows  formally  from 

j  9 

Corollary  1  to  Theorem  8  by  the  substitution  of  ip  for  t  where 
2 

i  = -1.  In  particular,  since 

i 

cosh  x  =  cos  ix,  sinh  x  -  -i  sin  ix 

we  have 

‘Vlp)  =  ^n(p^  4*n  (ip )  =  4>n(p) 

formally.  Precisely  the  same  type  of  proofs  used  for  the  sin  xt  and 
oos  xt  cases  establishes  the  exactness  of  (10.  6)  for  the  functions 
(5.  3),  where  here  we  need  the  sinh  and  cosh  functions  instead  of 
the  sine  and  cosine  functions,  respectively.  Where  before  (see  (8.  26) 
and  (8.  28))  we  had 

IHJ  s  | eitv |  +  K2Pj  V  jxj  >  1,  constants 

now  we  obtain 

|K  |  §  k  e  vp  +  K,  p.  V  p.  >  1,  K  ,  K  constants 

v  3  41  1*3'  4 

so  that  (24)  is  satisfied. 
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IV.  EXPRESSIONS  FOR  THE  ERROR 

11.  An  explicit  expression  for  the  remainder  Rf.  In  the  intro¬ 
duction  we  mentioned  that  there  was  still  a  third  approach  to  our  par¬ 
ticular  q.  f.  for  the  case  of  odd-degree  splines.  This  way  lay  through 
the  use  of  the  so-called  Rodrigues  function  H(x)  of  the  Peano  kernel 
of  the  q.  f.  [lOj.  We  consider  an  interval  of  integration  [a,  b]  and 

assume  n  interior  nodes  x  such  that 

v 

a  <  x,  <  . . .  <  x  <  b. 

1  n 

/  m  \ 

Let  w(x)  e  I.(a,  b)  be  a  given  weight  function  and  let  w'  *(x),  m 
a  fixed  integer  =1,  denote  any  2m- fold  integral  of  w(x).  Suppose 
I  U  I'  is  defined,  not  as  in  (5,  15)  and  (5.  16),  but  as  follows;  Let 
I  be  a  subset  of  {0,  1, . . . ,  m- 1),  I  =  {0,  1, . . . ,  m-1  }  -  I  and 
I‘  =  {2m-l-i  :  i«  1°}.  Define  J  U  J*  similarly.  Schoenberg  in  [  10, 
§?)  discusses  so-called  complete  quadrature  formulae  of  the  form 

(1!.  1)  /  w(x)f(x)dx  =  £cf(x|+  y  A(f(1)(a)+  £  Bf(J)(b)  +  Rf 

a  ►.!  v  v  faiur  1  j.jur  1 

where 

(15.2)  Rf  =  /  H(x)f(2m)(x)dx, 

a 

Under  suitable  conditions  on  the  sets  I  and  J,  q.  f.  of  the  form 
(11.  i)  and  (11.2)  exist  for  any  choice  of  weight  function  w(x).  In 
this  event  the  H(x)  of  (11.  2)  is  a  unique  monospline  of  the  form 
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01.3)  H(x)  =  w(“2m)(x)-  S_  (x) 

2  m,  n 

where  (x)  is  a  spline  function  of  degree  2m  -  1  with  the 

2  m,  n 

simple  knots  x^, . . . ,  x^,  satisfying 


o 

ii 

•n 

if  i  «  I  U  I' 

01.4) 

E(xJ  =  0 

v  =  1,  2, . . . ,  n 

H(j)(b)  =  0 

if  j  €  JU  J*. 

The  coefficients 

C  ,  A.  and  B. 
v  i  J 

are  given  by 

C  =  H(Zm"  \  -0)  -  H(2m“  %  +0)  v  =  l,  2, . . . 

01.5)  Aj  =  -(-l)1  H(2m"1'1)(a)  if  ic  III  V 

8^  =  {-i)J  H<2m~H)(b)  if  j  c  JU  J' 


y 


n 


We  want  a  related  expression  for  the  interval  [0,  ao)  and  the  fol¬ 
lowing  choices  of  weight  functions  w^(x)  and  set  I  U  I'  . 


(11.6)  wt(x)  =  c:OS  xt>  I  U  I*  =  {0,  1,  3,  5, . . . ,  2m-3  } 

(11.7)  wfc(x)  =  sin  xt,  IU  V  =  {0,  2,  4, ,  2m~4,  2m-2}. 


To  obtain  an  approximation  like  (11.  3)  for  H(x),  we  first  consider 

spline  interpolants  to  wt(x)  on  the  whole  line  R.  Because 

I  A2tn”  1  ,  ^  I  ,2m-  1 

sup  I A  wt(v)|  S  2  <«, 

v 

by  Theorems  1  and  2  of  [9,  p.  1 69 ]  we  know  that  there  exists  a  unique 

cardinal  spline  function  St(x)  satisfying 

(11,8)  SJv)  =  wjv)  for  all  integers  v 
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and 


(11.9)  S  (x)  <  S  ft  L2m_1(R). 

t  2m  oo 

On  the  other  hand,  since  the  sequence  {wt(v)  }  is  bounded,  by 
Theorem  1  of  [11]  there  also  exists  a  unique  cardinal  spline  function 
§t<x)  satisfying 

(11.10)  St(v)  =  wt(v)  for  all  v 
and 

(11.11)  §  (x)«  S  fl  L  (R). 

t  2m  oo 

We  want  to  show  that  St(x)  is  the  same  as  SJx).  From  the  nature  of 
the  data  (11.  6),  (11.7),  we  know  that 


(11.  12)  St(x)  and  §t(x)  are  even  or  odd  as  the  sequence  {wjv) }  is 
Let  S(x)  -  St(x)  -  §t(x),  so  that  (11.  8)  and  (11. 10)  imply  that 


(11.  13) 


S(”)  =  0 


for  all  integers  v. 


Evidently  (11.9)  and  (11.  11)  require  S(x)  t  S_  so  that  we  also  have 

2m 

S(x)  c  .  We  wish  to  show  now  that  S(x)  can  grow  by  at  most 
some  power  of  x,  so  by  (11.  12)  it  is  sufficient  to  consider  S(x)  for 
x  =  0.  In  particuU  t,  we  can  write  S  (x)  in  the  form 

......  m  .  i  .2  m- 1  2  m- 1 

(11.14)  S  (x)  =  a  +  a  x  +  . , .  +  — - —  x 

i  u  i  (2m- 1)! 


+  (2m- 1)!  ECM 

v  =  l 


2m- 1 


whore  ‘ho  coefficients  are  to  be  determined.  Taking  2m- 1  derivatives 
in  (9.  14)  gives  the  relations 
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St' "’1)(x>  *  a2m-l  +  %  Cv(x"’,)t 


S<2m-1|(x,  . 


a,  . ,  +  Y  C  if  x  is  an  integer 
2m+l  LJ,  v 
v=l 


i-  . ,  +  C  if  not. 
2m+l  -  v 
v  =1 

.  i _ i _ i  .t, e(2m~l). 


If  the  sequence  {C^}  is  not  bounded,  then  7(x)  /  L^fF) 

which  contradicts  (11.9).  Hence  we  have  the  following  necessary 


condition 


(11.15) 


s|2m_1)(x)  c  L^fF)  implies  {C^  >  is  bounded. 


Suppose  I  CM  <  K  for  v  5  ij  then  by  (11. 14)  we  have 


(11.16)  |St(x)|  ^  ^ao  +  aix+ ’•*  +(^T^!x2m~1| 


Si  .2m- 1 

(X**v)+ 


(2m- 1)! 


But  the  sum  in  (11. 16)  is  bounded  by  x 


2  m- 1 


so  that  we  get 


St(x)  =  0(x2m~S  as  x -*•  oo  . 

This,  with  (11. 11),  (11.  12)  and  the  definition  of  S(x),  implies  that 


(11.17) 


S(x)=0(|x|2m  l)  as  x  -  +  oo  . 


A  special  case  of  Lemma  2  of  Schoenberg's  [11,  p.  12]  states  that  if 


S(x)  «  s!,  and  satisfies  (11. 17),  then  S(x)  is  identically  zero.  Thus 
<.m 

§t(x)  =  S^fx)  and  we  have 


-78- 


#1183 


? 


£ 


r 


(11.  18)  s  (x)  €  S,  n  L  (R) 

t  2m  oo  ' 

In  terms  of  St(x),  we  define 

(11.19)  Ht(x)  ="^(wt(x)-St(x)). 

From  the  form  of  wfc(x)  in  (11.6)  and  (11.7)  and  by  (11.9)  and  (11.18), 
we  obtain 

(11.20)  H  (x)  £  L  (R)  n  L2m~1(R). 

I  00  00 

We  also  seek  the  analog  of  (11.  4)  for  our  half-line  case.  Indeed, 
since  S^(x)  c  C  and  (11. 12)  holds,  we  have 

(11.21)  S^O)  =0  if  i  €  I  U  I' 

which  implies,  by  (9. 19),  that 

dl.22)  H^fO)  =  0  if  i  c  I  U  I'. 

Evidently,  (11.  19)  also  enables  us  to  write 

(11.23)  H^(v)  =  0  for  all  v. 

We  can  state  the  following 

Ifreorem  10.  _ Suppose  HJx)  is  given  bv  (11. 19)T  and  wt(x)  and 
1  U  I»  are  as  In  (11.  6)  and  (li.  7).  If_ 

(11.24)  «x)  .  C2m(K+)  n  L2m(K+)  n  L,(R+) 
then 

00  00  , 

(11.25)  /  f(x)wt(x)dx  =  fj  H^^ffv)  +  Tj  B  f'*\o)  +  Rf 

0  v=0  V>1  U I UI'  1 

where 
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(11.26) 


Rf  =  /  ^(x)  f(2m)(x)dx  . 

0 

Proof.  Since  (11.20)  and  (11.24)  hold,  we  write,  letting 
H(x)  =  Ht(x), 


00 


(11.27)  /  H  (x)f(2m,(x)dx  =  lim  H  (x)f'“,,,,(x)dx. 

0  b  oo  0 

By  successive  integrations  by  parts,  we  find 


e(2m). 


(11.28) 


/  H(x)f(Zm,(x)dx  =  [Hf(2m-1)-H‘f(2m“2)  +  ... 
0 

+  (-l)2m“2H(2rn-2i,]|Cb,+  (-I)21"1”*!  /b3H(2m“l) 


(x)f’  (x)dx  . 


But  H^2m  1  \x)  is  a  step  function,  so  we  split  up  the  interval  of 


integration,  and  find 


(11.29)  -  /  3H(2m“1,(x)f(x)dx  =-  /  H(2m“1)(x)f 


(x)dx 


m- 1) 


v=0  v 
v+1-0  v4l 


v=0 


(x)f(x)  [  -/  f(x)H  2m  (x)dx] . 

v+0  v 

r(2m). 


We  note  from  (11. 19)  that  we  can  substitute  wjx)  for  H  (x)  in 
(11.  29),  so  that  after  summing  and  rearranging,  (11.  29)  becomes 

(11.  30)  -  {-H(2m'1)(0+0)f(0)  +  ^  £  [H(2m-1)(v-0)  -  H(2m_1)(v+0)]f(v) 

v=l 

+  H(2m_  1  )([b]-0)  f([b]} }  +  }  «x)w(x)dx. 

0 


If  we  substitute  (11.  30)  for  (1 1.  29)  in  (11.  28)  we  find  that  (11.  28) 
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upon  evaluation  of  the  term  in  brackets,  becomes 


(11.31) 


/  H(x)f(2m,(x)dx  =  J  (- 1  )vHCv)([b])f <2m“ J” v) 


«[bj) 


v=0 


-  Y,  (“l)VH{v)(0)f(2rn’1*‘v) 
v=0 


(0) 


+  (-l)f H(2m“ 1  >(<H-0)f (0)  +  ^  [H(2m“1)(v-0)-H(2m"1)(v+0)]f( 


v) 


v=l 


/h) 

f(x)w  (x)dx  . 

a  • 


By  (11.  31),  then,  (11. 20)  and  (11.  24)  give  us,  on  letting  b  -*«,  that 


00 


2m~2 


(11.32)  /  H(x)f(2m)(x)dx  =  -  fj  (-l)VH(v)(0)f(2m"1"v)(0) 

0  v=0 


oo 


+  H(2m-l)(0+0)f((J)-  £  (H(2m”1)(v-0)-H(2m‘,1)(v+0)]f(v) 

v=l 


00 


+  /  f(x)wt(x)dx  . 


If  we  use  (11.  22)  in  (11.  32)  and  solve  for  /  f(x)w^(x)dx,  we  obtain 


00 


(11.33)  /  f(x)w  (x)dx  =  Yj  (~l)VH(v)(0)f(2m"1~v)(0) 

0  v/lUl' 


00 


-  H^2m~1)(0+0)f(0)  +  Y  [H(Znl"1)(v-0)-H(Zm"1)(v+0)Jf(v) 

v=l 


00 


+  /  H(x)f  2m  (x)dx  . 
0 
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Upon  interchanging  the  order  of  summation  in  the  first  sum  on  the  right 
of  (11.  33),  we  get 


00 


(11. 34)  /  f(x)w  (x)dx  =  -  H 

0 


(2m“1)«H0)f(0)  -  Z  (-l)1H<2m"1"i)(0)f(i,(0) 


ielUP 


oo 


00 


+  l  [H(2mrl)(v-0)-H(2m"1)(v+0)3f(v)  +  /  H(x)f(Zm,(x)dx. 
v=l  0 


If  we  define 


h!2"0  =  -  H(2m“1}(0+0) 

L 

(11.35)  H*2"0  =  H(2m_1)(v-0)  -  H^^fv+O)  v  =  1,2,... 

w,  t 


B(2m)  =  -(-l)1H<2nh’1"1,(0) 
*> * 


i  «  I  U  I' 

we  obtain  the  desired  form  for  our  q.  f.  (11.  25)  and  (11.  26)  and  so 
establish  the  theorem. 


1 2.  The  remainder  Rf  for  the  cosine  transform  (5,  51).  We 

specialize  to  the  case  (11.  6)  and  now  want  to  establish  that  the  q.  f. 

in  Theorem  10  is  the  same  as  the  q.  f.  in  Theoi  -  7.  We  do  this  by 

examining  the  form  that  the  function  Ht(x)  as  defined  by  (11. 19) 

must  take.  Enforcing  the  requirements  (11.  22)  and  (11.  23)  will  lead 

us  by  a  generating  function  approach  to  tne  same  coefficients 
(2m)  D(2m) 


H 


7,  ,  of  Theorem  7. 

',  t  2J-1,  t 


We  attempt  to  find  an  expression  for  the  cardinal  spline  function 
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St(x)  satisfying  (11. 8)  and  (11.9) .  Since  we  need  only  consider  the 
half- line  F+,  St(x)  can  be  written  in  the  form  (11. 14)  where  the 
coefficients  a^,  C,  are  to  be  determined.  We  want  to  express  the 

in  terms  of  H^2™^  =  H  and  =  B.  .  From  the  definition 

of  Ht(x)  (11. 19)  or  from 


(12.1) 


(-l)mt2mHt(x)  =  cosxt-St(x) 


we  find  by  differentiating  2m- 1  times  and  employing  (11.  29)  that 
(12.2)  Cy  =  (-l)mtZmHv,  a2m_1  =  (**l)mt^mH0  . 


If  we  differentiate  (12. 1)  (2rn-  1-i)  times  and  use  (11.  29),  we  get 

2m-l-i 


(-l)mt2mB  =  (-1)  2  t2”1"1’1  _  a 

1  2m- 1-i 


i  £  I  U  I' 


or  solving  for  a,  .  .  that 

2m'1'1  ^-i-i 

(12.3)  a2m-l-i  r  (‘U  2  (-l)mt2mB1  1,1  UP. 


Enforcing  (11.  22)  in  (12.  1)  gives  us 


(12.4) 


at  =  0  if  i  c  I  U  I’ 


so  that  by  employing  (12.  2),  (12.  3)  and  (12.  4)  we  may  rewrite  (11.  14) 


[-t2-(-l)mt2mBZn1_,]  2  _  [t’-(-l,VmB2m.5)  4 


.4  ,  ..mini. 


(12.  5)  SJx)  =  1  4 


...  +  [(-I)m“1t7m-2-(-i)mt2mB13 

i  (2m- 2)! 
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a  *~»**w*<^»  j~*»  «  -  ,  ,  ««  —  ^  r 


[mA2i\1  n)Vm  f  to-i 

(2m- 1)!  +  (2m-l)!  ^  v  4 


By  virtue  of  (12.  2)  and  (11. 15)  we  have  the  following  necessary  condi¬ 
tion: 

(12.6)  sj2m-1)(x)c  L  (R)  implies  {Hv }  is  bounded. 

If  we  solve  (12. 1)  for  cos  xt  and  then  require  (11.23)  for  k  a  positive 
integer,  we  obtain  the  sequence  of  relations 


(12.7)  cos  kt  =  1  +  Yj  [(-l)lt2j-(-i)mt2mB2rn^i_2j] 

1 —1 


j CL 

(2j)! 


,  ..m2m  oo  -  . 

+  izlJL_t -  y  H  2m- 1 

(2m- 1)!  A  v'  4 


(k  =  1,2,...). 


The  form  of  these  relations  suggests  the  use  of  generating  func¬ 
tions  for  the  determination  of  the  coefficients.  The  righthand  side  of 

k-1 

(12.  7)  in  view  of  Lemma  3  is  equal  to  the  coefficient  of  x  in 


j  m-1  [(-l)Jt2J-(-l)mt2mB  }  P  (x) 

(  *  1-x  +  A  (2j)!  .21+1 


(1”X) 


,  ..m^m  oo  P„  _(x) 

+i7^TTTr«  l 

v^0  v  (,.x,2m 


The  lefthand  side  of  (12.  7)  is  the  coefficient  of  x  in  what  we  called 
t(x),  defined  in  (4,39).  Equating  t(x)  and  the  expression  in  (12.8) 
and  then  solving  for  £  gives 
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(12.9)  l 

v-0  V  t2m-2w  \  rm 


m“1  hl)it2i  r  ,  .2m- 1-21. 

4 — 1 -  P-.  .(x)(l-x)  '] 


in 

-s 


i=l 


(2J)!  2j~f 


m-l  B„  ,  ..  _ 

J-!  (2j)!  2J- 1 '  ni  x) 


A  change  In  the  order  of  summation  in  the  last  sum  on  the  right  side 
of  (12.9)  leads  to  the  final  relation 


00 


v  ..  v  (2m-l)!  (-1)  r.  .,2m  .  .  ,  ..2m-l 

(12.10)  l  H  x  =  -J - 7UT<  t(x)  +  (x-1) 


v=0 


P2m-2(X)1  t2m 


.  1  t2^  r  .  ..  2m-l-2J. 

+  j-l  1  (2))!  P2j-1  X  X"*1  ^ 


m-l 


?.  B2j-1 


j  =  l 


(X-1?2i"lp2m-l-21(X)> 
(2m-2j)> 


This  is  precisely  the  same  relation  as  (5. 1)!  The  analysis  in  §5 
led  to  a  unique  choice  of  the  sequence  {H^  }  under  the  stipulation  that 
this  sequence  be  bounded.  By  (12.  6)  and  the  existence  and  unicity  of 
an  St(x)  satisfying  (1 1.  8)  and  (l  l.  9),  we  conclude  that  the  and 
Bj  as  determined  in  §5  are  the  required  coefficients  for  S^(x)  as  given 
in  (12.5).  So  this  approach  through  the  use  of  the  function  H^lx)  leads 
to  precisely  the  same  q.  f.  as  that  of  Theorem  7,  and  in  particular  leads 
to  an  expression  for  the  remainder  Rf  in  Theorem  "  We  have  thereby 


established 

#1183 


-85- 


•> 


Theorem  11.  The  remainder  Rf  in  the  q.  f.  (5.  51)  of  Theorem  7 
under  the  stipulation  (5.  52)  for 

(12.  11)  f(x)  €  C2m(R+)  n  I2m(R+)  n  Lj(R+) 

is  given  by 

(12,12)  Rf  =  /  H.(x)f^2m)(x)dx  =-M^r  /  [cos  xt  -  S  (x)]f*2in\x)dx 
0  t  o 

where  St(x)  is  the  unique,  bounded  (2m- l)st  degree  cardinal  spline 
interpolating  cos  xt  at  the  integers. 


13.  A  bound  on  the  remainder  Rf  of  (5.  51).  We  now  examine  the 

problem  of  expressing  the  above  cosine  transform  (5.  51)  in  steps  of 

length  h.  Let  f(x)  c  S_  n  L.(R+)  so  that  by  Theorem  7  and  (5. 42) 

2m  i 

we  can  write 


00  ij>  (t)  f  oo 

(13.1)  /  f (x)co s  xt  dx  ■■  Mf(0)  +  Y.  f(v)cos  vt 

0  *2mW  |2  v=0  1 

"K1  l  ill  -M(t>  U 

+  i  \  2:~  j  ) f(zj  t>(0). 


+2m(t> 


J=1  t 

Let  now  F(x)  be  a  (2m-l)st  degree  spline  function  in  (0,  oo)  having 
its  knots  in  x  =  h,  2h,  3h, . . .  where  h  >  0.  We  want  to  express  the 
cosine  transform  in  terms  of  the  values 

(13.2)  F'(0),  F'"(0),...,F(2m'3)(0),F(0),F(h),F(2h),...  . 

If  we  let 


(13.3) 


f(x)  =  F(xh) 


(0  I  X  <  oo) 
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we  see  that  f{x)  is  a  semi-cardinal  (2m-l)st  degree  spline  function 

i 

for  which  the  data  are 

(13.4)  f«(0)  =  hF'(0),...,f(2m“3)(0)  =  h2m"V2m~3)(0), 

f(0)  =  F(0),  f(l)  =  F (h),  ...  . 

From  (13.  1)  we  therefore  obtain  the  relation 


80  tp  (t)  f.  00  > 

(13.5)  /  F(xh)cosxtdx  =  •  J“F(0)+  T.  F(vh)cos  vt  [ 

0  ^2  nr  J  2  v=l 


+  £  -M^-H  - -2i 2,?,'21  1  h2i  1  F(2ill(0). 

J=1  t2)  *2m(t) 

Replacing  in  the  integral  x  by  x/b  and  replacing  afterwards  in  this 
identity  t  by  th,  we 

^ro(th> 


00 

(13.6)  /  F(x)cos  xt  dx 


*2m{th> 


hi  “ F(0)  +  F(vh)cos  vth] 


v=0 


>i(thH2m-2j(th>  ,„(2J-1) 


j  =  l  t 


2J 


<*>2m(th) 


(0). 


Suppose 

(13.7) 


.  _,2m  _  2  m _  _  (T,+. 

F(x)  (  C  n  Lj  (R  )  n  Lj(R  ) 


then  (13.  6)  is  no  longer  an  identity.  However,  the  righthand  side 
will  give  us  the  desired  approximation  to  the  cosine  transform  of  F(x) 
for  reasonably  small  h.  We  now  want  to  see  how  good  an  approxima¬ 
tion  this  is. 
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(13.  5)  is  still  valid  if  we  use  (12.  15)  and  (13.  3)  and  add  the  term 

i  m  00 

(13.  8)  .  h2m  /  [cos  xt  -  S  (x)}F{2m)(xh)dx 

t  0 

to  the  right  side  of  (13.  5).  If  we  denote  (13. 8)  by  R  and  again  first 
replace  x  by  x/h  and  then  t  by  th,  we  obtain  the  expression  (13. 6) 
with  the  added  term 

(13.9)  R  [cos  xt  -  Sth(xA)]F^m,(x)dx 

on  the  right  side.  Here  S^(x/h)  is  the  unique,  bounded  (2m-l)st 
degree  spline  interpolating  cos  xt  at  vh  for  all  integers  v.  We 
get  a  bound  on  R  by  bounding  what  we  shall  call 

(13.  10)  M(t,  h)  =  max  .  [cos  xt  -  Sth(x/h)]  * 

x  =  0  t^m 

♦ 

X  2it  z 

Let  z  so  that  SJL.(— r~)  is  the  unique,  bounded  (2m- l)st 

2w  th  h 

h  ?h 

degree  spline  interpolating  cos  2*zt  for  z  =  0,  +  — ,  +  ~  ,  . .. 

We  also  let  h  be  of  the  form 

2 

(13. 11)  h  =  — -  for  n  a  natural  number 

n 

so  that  the  spline  agrees  with  cos  2iTZt  for  z  =  0,  +—  ,  +  —  ,  ... 

7  —  n  ~  n 

If  we  consider  t  =  1,  2, . . . ,  n-1,  then  cos  2?tzt  is  periodic  on  the 
interval  [0,1]  and  so  is  S^(2TTz/h),  by  Theorem  6  of  [11].  We  require 
a  special  case  of  Lemma  6.  3  of  Golomb' s  [3]  which  we  state  as 


-88- 


#1183 


wm 


w;v 


§ 


Lemma  6  (Golomb).  Let  bg(u)  denote  the  unique  bounded 
(2m-l) st  degree  spline  that  Interpolates  the  function  e2*lsu  ^ 
k/n  (k  =  0,±l,  +  2,... ).  Then 

(13  12)  |  2irisu  ,  .  j  2m  2m  -2m 

U3.12)  |e  bg(u)  1=4*2  s  n  s  =0,  +1, . . . ,  +(n-l). 

In  [3,  p.  13]  Golomb  remarks  that  Re  b  (u)  and  Im  b  (u)  are  the 

s  8 

corresponding  spline  Interpolants  to  cos  2ttsu  and  sin  2ttsu, 
respectively,  so  that  we  also  have 

03. 13)  |cos  2ttSU  -  Re  bs{u)!  *  4-  22mS2  V2m  Su-o, +1, . . .  ,±(n-l) 

where  Re  bg(u)  is  the  unique,  bounded  (2m-l)st  degree  spline  that 
interpolates  cos  2ttsu  at  the  points  k/n  (k  =  0,  jfl,  +2,  . . .  ). 

We  can  therefore  employ  (13. 13)  to  get  the  bound 

(13.14)  max  |cos  2T?tz  -  St,  (2wz/h)[  S  4-  22rnt2mn”2m 
z  =0  th 

^  ~  2  j  •  •  *  |  1  • 

Combining  (13.14)  and  (13.11)  with  (13,10)  and  the  definition  that 
z  =  x/2ir,  we  obtain 


2 

(13.15)  M(t,h)  =  max  i~  [cos  ZTrtz  -  S,.  (2nz/h)l|  ?4~ 

z£0  t  th 


2  m 
2m 


t  =  1,  2, . . . ,  n-1 . 

We  now  consider  t  in  the  form 

(13#  r6}  t-  =  p/q  where  p,  q  t  Z  and  p  =  1,  2, . . . ,  n-1 . 

Now  let  w  =  z/q  and  we  get  from  the  equality  in  (13. 15)  that 


#1183 


-89- 


(13.171  M(t,h)  =  max  l-L-  (cos  2„pw  -  fi.k(^)j| 

VbO  t2m  th  h 

_  ,2trqw. 

where  sth  (*“£[—)  interpolates  cos  2*pw  for  w  =  0,  +  h/2«q, 

+  2h/2trq, ...  If  we  let  h  be  of  the  form 

(13.18)  h  =  2nq/h  for  n  a  natural  number 

and  consider  p  =  1,  2, . . . ,  n-1,  then  Sth  (-^p )  is  periodic  for 
0  =  w  =  1.  By  (13. 13)  then  we  have 


(13. 19)  max  jcos  2trpw  -  S^(^^)|  5  4-  22mp2V2m 

w  >  o  th  h 

(p  =  1,  ...,n-l). 

Substituting  (13. 19)  in  (13. 1?)  in  view  of  (13. 16)  and  (13. 18)  we  find 

,2m 

(13.20)  M(t,h)  S4^ 

Suppose  now  we  fix  h  of  the  form 


(13,  21)  h  =  ~  for  N  a  natural  number. 

We  let  q  be  any  positive  integer  and  choose  n  such  that 


(13.  22) 


3.  1 

n  ~  N  * 


Then  (13,  22),  (13. 18)  and  (13.  20)  enable  us  to  show  that 


(13.23) 


,2m 

M(t,  h)  =  4  - — 
2m 

IT 


<t=M. 
q  q- 


By  (13. 19),  (13.10),  (13.21)  and  (13.23), 
(13.24)  I R |  §  M(t,h)!!F(2m)|L  5 


then,  we  have  the  bound 
,2m 


2m 


||  p(2m)  jj 


(t.i  1  1H  1  V 

q  q  n  q 
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and  the  following 


Theorem  12,  Suppose  ft  C2m(P+)  n  L.(F+)  n  L?m(R+)  and  — 

1  l  ■  n 

is  a  natural  number.  Then  we  can  bound  f  Rf  ]  as  given  in  the  q.  f. 

r°  t  »  > 

(13. 25)  /  f(x)cos  xt  dt  =  - — h  -  f(0)  +  £  f(vh)cos  vtlA 

0  v2m'  v=l  J 


vT1  (-nJ  <Wth^2m-2^th*  m  n 

♦  l  -L£n  -  -iJ  -.  -  zr121  ]f(2HI«»  * 

)=1  r*  *2m<U” 


Rf 


by 


(13.26)  I  Rf  I  *  4(^)2m|Jf^2m^||  +  for  all  rational  t  in  (O,-^). 

»  ij.iK  /  n 


W.  ,^g^!g|e^^2rem!^_of_J2.  We  also  have  an  analogous  theorem 
to  Theorem  1 1  for  sin  xt,  which  we  shall  only  state  as 

Theorem  13.  The  remainder  Rf  in  the  q.  f. 

»j<  (t)  foo  > 

(14.1)  /  f(x)sin  xt  dx  =  .  J  Y  f(v)slnvt' 

0  ^2m  *  \v=l 

♦  21  ■%£  (1  -  ^]f(2»(0)  +  Rf 

)=0  t2,+1  ♦zm'1'  * 

fgr.  f  c  C2m(F+)  n  lJ2it^(F+)  n  L  (F+)  is  given  by 
t  n®  00 

(14.2)  Rf  =  /  [sin  xt  -  S  <x)Jf(2m,(x)dx 

o  8 


where  g(x)  is  the  unique,  bounded  (2m-l)st  degree  cardinal  spline 
interpolating  sin  xt  at  the  Integers. 
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Let  S{  c(x)  denote  the  spline  Sfc(x)  of  Theorem  11  interpolating  cos  xt. 
Toward  a  proof  of  Theorem  4,  let  us  assume  that 


(14. 3) 

Then  we  can  write 


f «  c2m  n  L2m(F)  n  Lj(R). 


If'  W 

(14.4)  J  f(x)e  dx  =  /  f(x)cos  xt  dx  +  i  /  f(x)sinxtdx 


*  w 
=  J  [f(x)  +  f  (-x)Jcos  xt  dx  +  if  [f(x)  -  f(-x)]sin  xt  dx. 
-co  0 

Applying  (13. 25)  for  h  =  l  and  (14.  1)  to  f(x)  +  f(-x)  and  f(x)  -  f(-x), 
respectively,  we  find  that  (14.4)  becomes 

<H.s,  K!' 

*2mW  vr-oo  j=l  (tZ))  *2m(t) 

•  [f(2,“1)(0+)-  f^"  1>(0" )J  +  R  f 


{»|»_  (t)  »  -00 

TfTil  2jf(v)sinvt+  Z  f(v)sinvt] 

"2  m'  v  =  l  v=- 1 

*"%  ^][f,2))(0+)-f<2"(0-)]  rs4 

«  t  “2m' 


J=0  f 


where 


co  m 

(14.6)  Rcf  =  f  “  [cosxt-  S^c(x)If(2m(x)  +  f(2m,(-x)]dx 


oo  ,  ..m 


f  =  /  [sin  xt  -  S  (x)Xf(2m)(x)  -  f<2m)(-x)Jdx. 

o  t  8 


We've  taken  ^(t)  -  1  for  convenience  of  notation.  By  (14.  3)  the 
derivative  terms  in  (14.  5)  are  zero  so  that  substituting  (14.  5)  in  (14.4) 


gives 


(14.7) 


r° ..  .  ixt .  ^2m^  V  ...  ivt  , 

J  f(x)e  dx  =  ~ — rr  2j  *<v)®  +  Rcf  + 

-co  ^2m'  v=-oo 


From  (14. 6) 


«?  .  ..m 


R  f  =  /  [cos  xt  -  S.  rj(x)Jf(2m)(x)dx 

c  '  t2m  t,  c 

+  /  ^-[cosf-ty)  -  S  (-y)]f(2m)(y)(-dy) 

0  t2m  l» c 

where  we've  used  y  =  -x  to  get  the  second  integral.  But 

st  =  S.  (y)  so  that 
t,  c  t,c 

(14. 8)  Rcf  =  /  [cos  xt  -  St  c(x)]f*2m)(x)dx. 

-00  t  * 

We  can  argue  the  same  way,  using  the  fact  that  S  (-y)  =  -  S 

s  t, 


and  obtain 


(14.9)  Rgf  =  /  [sinxt-  St  s(x)]f(Zm)(x)dx. 

-00  t  * 

By  (14. 7),  (14.  8)  and  (14.9),  the  remainder  Rf  of  Theorem  4  therefore 


has  the  form 


°c  .  ..m 


(14.  10)  Rf  =  Rcf  +  iRgf  =  /  [elxt-(St  c(x)  +  iSt  s(x))]f(2m|x)dx. 

-00  .  t  *  * 

We  note  that  S  (x)  +  iS  (x)  is  the  unique,  bounded  (2m-l)st  degree 
t,  c  t,  s 

ixt 

cardinal  spline  that  interpolates  e  at  the  in*  gers.  By  precisely 
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the  nine  argument  that  led  to  the  bound  (13. 26)  of  Theorem  10  except 
that  now  we  use  (13.  12)  of  Lemma  6  where  before  we  used  (13.  13), 
we  reach  the  desired  bound  (2. 17)  of  Theorem  4. 

In  the  foregoing  proof  in  ( 14.  10)  we  determined  the  explicit  form 
of  Rf,  a  result  we  state  as  a 

Corollary.  For  f(x)  c  C2m  n  L2m(R)  n  Lj(R),  the  remainder  in 
the  q.  f. 

(14. 11)  /"  Hx)ebttdx  =  £  «v)e‘vt  tRf 

-»  *2m(t'  v=-« 

is  given  bv 

<14. 12)  Rf  =  /”  [e1*1  -  b  (x)J  f(2m)(x)dx 

t  -00  1 

where  b^fx)  Is  the  unique,  bounded  cardinal  spline  interoolant  of  the 
function  e***  at  the  integers. 


15.  An  explicit  expression  for  the  remainder  Rf  in  Theorem  3.  The 
same  type  of  approach  we  took  in  §11, 12  will  lead  us  to  an  explicit 
form  for  the  remainder  Rf  in  Theorem  3.  For  our  case  of  the  weight 
functicn  e  and  I  U  I'  =  {0, 1, . . . ,  m- 1 }  we  require  the  unique 
semi- cardinal  spline  S  (x)  <  S,  n  L_(F+)  satisfying  the  conditions 

p  cm  l 


(15.  1) 
and 


S  (v)  =  e 
P 


-  «TVP 


(v  =  1,  2, ... ) 


Here  (15.2)  reflects  the  fact  that  the  first  m-  j  derivatives  of  e~Xp 
and  Sp(x)  agree.  In  terms,  of  SJx)  we  can  state  the  following 


Theorem  14.  Let.  f(x)  *  C2m  and  f(i)(x)  t  F  for  some  v  . 

Yj  -  T 

J  =  0,  1, . . . ,  2m- 1.  Then  the  remainder  Rf  of  Theorem  3  may  be  ex¬ 
pressed  as 


(15.3) 

where 


oo 

Rf  =  /  H  (x)f,2m,(x)dx 


<l5-4>  Hp,x)  =  l6‘XP  - 

P  K 

In  order  to  obtain  the  appropriate  form  (23)  of  Theorem  3  we  want 
an  analog  of  Theorem  10  to  hold.  By  (15.  l)  and  (15.  2),  we  find 
(15.  5)  H*j)(0)  =0  if  j  €  I  U  I* 

and 


(15.  6)  Hp(v)  =0  (v  =  1,  2, . . , ) 

so  that  an  appropriate  analog  can  be  proved  in  the  same  way  that 

Theorem  10  was  proved  if  we  can  show 

(15.7)  lim  H*o)(b)f(P)(b)  =  0  if  aftps  2m-l. 

b -oo  p 

By  the  assumptions  of  this  theorem  and  (15. 4),  this  amounts  to  showing 
(k) 

(15.  8)  8^  (x)  for  k=0, 1, . . . ,  2m-l  is  of  exponential  decay 

as  x  —  oo. 

Just  as  in  §8,  we  shall  discuss  only  the  cubic  case  m  =  2$  the 

extension  to  higher  m  is  very  similar.  By  using  (8.  6)  we  can  write 

S  (x)  as 
P 

#1183 


-95- 


CO 

(15.9)  S  ix)  =  2  e"vpL(x)-  pA(x) 

"  v=0  v 

where  Mx)  and  A(x)  are  the  fundamental  functions  discussed  in  §8. 
By  using  (8.  15),  (8. 21)  and  (8.  23)  in  (15.9),  we  obtain 

(15. 10)  S  (x)  =  L(x)  +  £  e‘vpL(x-v)  +  J  e"vp(VJ  XV)S_(x)-  -£■  S,(x). 

p  v=l  v— l  1  V3  1 

But  for  the  special  case  m  =  2  of  Theorem  2  of  [  11,  p.  5  ],  we  have 

i  i  “Y*^ 

(15. 11)  |L(x)|  <  Ce  for  all  real  x, 

where  C  and  y*  are  positive  constants.  This,  in  conjunction  with 
(8.  13),  and  the  expression  <15.  10)  guarantees  that  S^(x)  decays 
exponentially  as  x-*  ce. 

Similarly,  by  differentiating  each  side  of  (8.  10)  and  (8. 12)  three 
times,  we  obtain 

(15.12)  L’"(x)  =  -6^(X-l)Xk  if  k  <  x  <  k+ 1,  k£l 

and 

05.13)  Sj'" (x)  =  -6(X-l)Xk  if  k  <  x  <  k+ 1 

so  that  S^"(x)  as  determined  from  (15.  10)  also  decays  exponentially 
as  x  •*  oo.  Therefore  (15.  8)  holds  and  we  obtain  the  desired  analog 
to  Theorem  10. 

Now  we  employ  the  same  procedure  used  in  §12  to  show  that  the 
q,  f.  we’ve  obtained  is  indeed  the  same  q.  f.  as  given  by  Theorem  3. 

The  approach  in  §12  leads  us  to  the  relation  (9. 13)  so  we  can  con¬ 
tinue  as  in  the  analysis  of  §9  to  finally  obtain  our  desired  q.  f.  and 
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in  this  way  a  proof  of  Theorem  14.  The  actual  procedure  is  too  similar 
to  repeat. 


16.  Computational  examples.  In  [2],  Einarsson  compares  several 
methods  for  computing  cosine  transforms  for  the  special  case  of 

»x 

f(x)  =  e  .  One  method  he  uses  and  the  reason  for  the  paper  is  based 
on  the  approximation  of  f(x)  by  its  cubic  spline  approximation.  This 
q.  f. ,  precisely  the  same  one  as  (13.  25)  for  m  =  2,  is 


00 


(16.  1)  /  f(x)oos  xt  dx  = 

0 


Vth1  r  i 

*  ♦4<th)h|l 


00 


“  f(0)  +  Yt  f(vh)cos  vth^ 
v  —  1 


,  ^(th) 


where  we've  used  <{>2(t)  =  1.  Einarsson' s  main  conclusion  is  that  this 
spline  q.  f.  is  superior  to  Filon's  formula,  a  q.  f.  based  on  approxi¬ 
mation  of  the  function  by  a  quadratic  in  each  double  interval  and  one 
of  the  most  used  formulae  for  the  calculation  of  Fourier  integrals. 

One  of  the  other  methods  Einarsson  uses  for  comparison  is  the  so- 
called  Filon- Trapezoidal  rule  found  in  Tuck  [  14  j,  which  for  the  inter¬ 
val  (0,  oo)  is  merely  (13.  25)  for  m=l,  that  is,  the  linear  spline  case. 

Einarsson' s  calculations  indicate  that  for  small  values  of  t,  the 
q.  f.  (16.  1)  gives  a  relative  t.-ror  that  is  four  times  less  than  the 
Filon  formula.  For  large  values  of  t,  the  relative  error  of  the  Filon 
formula  increases  rapidly,  while  the  spline  method  (16.  1,  gives  a 
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surprisingly  small  error  growth.  This  same  phenomenon  we  found  to 
be  the  case  for  the  following  q.  f. ,  obtained  fcom  (13. 25)  and  (8.  3), 
respectively,  for  m=3 


(16. 2)  f  f(x)cos  xt  dx  = 
0 

,  «Mth) 


XJi 

Vth' 


f(0)  +  2]  f(vh)cos  vthj 

V=1 


i  ,  t.(th)iMth) 

+  7 11  "  ]f,(0)  ♦  Rf, 


(16.3) 


00 

/  f(x)cos  xt  dx  = 
0 


*6<th) 

♦6'(th) 


h 


00  ' 
f(0)  +  ^  f(vh)cos  vth  , 

V  —  1 


+4(th) 

«i>6(th) 


,  f4(th)  4«.(th) 

]f'(0)  +^[1  8"'(0)+  HI 


two  q.  f.  corresponding  to  spline  approximation,  the  first  using 
I  U  I'  =  {0, 1,  3  }  and  the  second  using  I  U  I'  *  {fi^  i,  2  }. 

We  now  consider  the  absolute  error  and  concern  ourselves  with 
the  q.  f.  (16.  1),  (16.  2)  and  (16.  3)  and  the  bounds  we  obtained  for  the 
error  in  (16.  l)  and  (16. 2)  for  two  examples.  We  first  remark  that 
as  the  step  h  gets  small,  it  appears  that  $"'(0)  of  (16.  3)  approaches 
f'"(0)  of  (16.  2)  so  that  the  difference  in  these  approximations  be¬ 
comes  very  small.  An  instance  of  this  we  indicate  below. 

In  Figures  1-6  the  absolute  values  of  the  absolute  error  fo.-  the 
calculation  of  the  cosine  transform  with  (16.  l)  and  one  of  (16. 2)  or  (16.  3) 
is  given  as  a  function  of  t  for  the  stepsize  h  at  80  different  places 
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from  t  ~  .  25  to  i  =  50.  The  cubic  curve  is  represented  by  x’s  and 
the  quin?:-  curve  for  (16.  2)  or  (16.  3)  is  solid,  the  dips  are  at  points 

i 

where  tn_-  absolute  error  changes  sign.  We  also  point  out  that  each 
axis  i*  soiled  logarithmically  and  a  lower  bound  for  the  error  in  each 
graph  is  2  x  10  .  Along  the  vertical  axis  we  indicate  by  a  3  or  a  5 

■  i 

where  the  computed  error  bounds  of  (13.  26). fall. 

In  Figure  l,  we  consider 


w 

/  e“X  cos  xt  dx  = 


1  +t 


and  the  stepslze  h  w  »2*  Here  S"'(0)  =-.99994  versus 
f*"(0)  *  -1  and  since  the  corresponding  graph?  arising  from  (16.  2)! 
and  (16.  3)  were  virtually  indistinguishable  we  only  need  consider  one, 
(16.  2).  We  also  compute  by  (13. 26)  the  bounds  on  |  Rf  |  and  find  that 

(16. 4)  I  Rf !  3  £  6. 1  x  10  5,  |  Rf  |  5  m  2.  4  x  10  7  forall  rational 

t  in  (0,  32) 

where  the  subscript  indicates  the  case  (16. 1)  or  (16.  2),  respectively.. 
We  note  that  the  bound  for  I  Rf  I .  is  actually  less  than  the  computed 
transform  corresponding  to  the  cubic  case. 


In  Figures  2  and  3  we  consider 


(16.  5) 


00  . 

/-*- 
0  1  +  x 


,  .  w  -t 
”  cos  xt  dx  =  —  e 


2tr 

for  h  =  «  .  4  to  examine  the  difference  between  the1  q.  f.  (16.  2)  and 

(16.  3).  We've  plotted  the  cubic  case  (16. 1)  also  to  serve  as  a  reference. 
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ERROR 


ABSOLUTE  VALUE  OF  ABSOLUTE  ERROR 


Figure  2.  f(x)  =  — l—  ,  h  =  7*  . 

1  +  x2  16 


cubic 


ERROR 
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VALUES 


We  see  that  for  small  values  of  t  the  q.f.  (16.  3)  (Figure  2) 
gives  a  considerably  worse  approximation  than  (16. 2).  Here 
S'M(0)  »  29?  compared  with  f"'(0)  =  0.  From  (13. 26)  we  find 

(16.  6)  I  Rf  I  j  »  1.  0  x  10  2,  |  Rf  1 5  2  5.  86  x  10  3  for  all  rational  t 

in  (0,16),  h=-^. 

In  Figures  3,  4,  and  5  we  again  consider  (16.  5),  but  now  for 
2^  2  2 

h  and  ■gj  respectively  to  consider  the  q.  f.  (16. 1)  and 

(16.  2)  as  h  decreases.  From  (13.  26),  we  obtain 

(16.7)  |Rf|3S  6.214  x  10-4,  |  Rf  |  g  £  9. 16  x  lo“5 

for  all  rational  t  in  (0,  32),  h  =  ~ 

(16.8)  |Rf|  3  S  3.884  xio"5,  |Rf|  S  1.  43  x  lo“6 

for  all  rational  t  in  (0, 64),  h  =  ~~  . 

The  figures  indicate  that  we  do  seem  to  have  errors  of  order  h*  for 
(16. 1)  and  h^  for  (i6. 2). 

If  we  consider  Figures  5  and  6  where  we  calculate  (16.  5)  for 
^  =  64  *rcm  2)  and  (16.  3),  respectively,  we  can  see  how 

much  closer  these  qulntic  curves  are  than  they  were  in  Figures  2  and 
3,  From  (16.  3)  tor  h  =  •“  ,  we  find  that  S"’(0)  = -.  00674  compared 
to  i!M(0)  -  0  and  th<*  S’ M(0)  -  295  we  computed  above. 


ERROR 


ERROR 
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